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Graphene is a famous realization of elastic crystalline two-dimensional (2D) membrane. Thermal 
fluctuations of a 2D membrane tend to destroy the long-range order in the system. Such fluctua¬ 
tions are stabilized by strong anharmonicity effects, which preserve thermodynamic stability. The 
anharmonic effects demonstrate critical behaviour on scales larger than the Ginzburg scale. In par¬ 
ticular, clean suspended flake of graphene shows a power-law increase of the bending rigidity with 
the system size, >c oc , due to anharmonic interaction between in-plane and out-of-plane (flexural) 
phonon modes. We demonstrate that random fluctuations of membrane curvature caused by static 
disorder may change dramatically the scaling of the bending rigidity and lead to a non-monotonous 
dependence of on L. We derive coupled renormalization-group equations describing combined flow 
of and effective disorder strength b, find a critical curve b{>c) separating flat and crumpled phases, 
and explore the behavior of disorder in the flat phase. Deep in the flat phase, disorder decays in a 
power-law way at scales larger than the Ginzburg length which therefore sets a characteristic size 
for the ripples—static out-of-plane deformations observed experimentally in suspended graphene. 

We hnd that in the limit L ^ oo ripples are characterized by anomalous exponent 2r] in contrast to 
dynamical fluctuations governed by rj. For sufficiently strong disorder, there exists an intermediate 
range of spatial scales where ripples decay much slower, with exponent r;/4. In the near-critical 
regime, disorder first increase with L, then reaches a maximum and starts to decrease. In this case, 
the membrane shows fractal properties implying a multiple folding starting from a certain length 
scale Li and finally flattens at a much larger scale L 2 (which diverges at criticality). We conclude 
the paper by a comparison of our results with available experimental data on graphene ripples. 

PACS numbers: 72.80.Vp, 73.23.Ad, 73.63.Bd 


I. INTRODUCTION 


Graphene, a single monolayer of graphitehas at¬ 
tracted enormous interest in the last decade (for review, 
see Refs. [i®. From the fundamental point of view, this 
interest is largely motivated by quasirelativistic charac¬ 
ter of its spectrum: charge carriers in graphene are two- 
dimensional (2D) massless Dirac fermions. This leads 
to a variety of remarkable phenomena. In particular, 
graphene is a unique example of a system where essen¬ 
tially quantum phenomena such as the quantum Hall ef¬ 
fect can be observed up to the room temperatureIn 
view of applications, the technological breakthrough in 
fabrication of flat monolayer 2D crystals offers new op¬ 
portunities in the future nanoelectronics. Remarkably, 
suspended graphene demonstrates extremely high room- 
temperature mobilityi^^— as high as 1.2-10® cm^/Vs and 
therefore is considered as a most perspective candidate 
for the carbon-based nanoelectronics. 

Elastic properties of graphene are also quite amazing. 
Free-standing graphene is an outstanding example of an 
elastic crystalline two-dimensional (2D) membrane with 
a high bending rigidity x ~ 1 eV. The most impor¬ 
tant feature distinguishing such a membrane from con¬ 
ventional 2D semiconductor systems is the existence of 
specific type of the out-of-plane phonon modes, so called 
flexural phonons In the harmonic approximation the 


energy of out-of-plane deformation reads 



dx 


ph^ -b K{Ahy 


( 1 ) 


where h(x) is out-of-plane distortion and p is the mass 
density per unit square. From Eq. o we find frequency 
of the flexural phonons 

uJc^ = Dq^, ( 2 ) 


where D = \Jxjp. Hence, in contrast to in-plane acoustic 
phonons, whose frequency scales as q, the flexural mode 
is very soft and, consequently, the out-of-plane thermal 
fluctuations of h{r) are very large. As a consequence, 
flexural phonons serve as a very effective scattering mech¬ 
anism for electrons (for discussion of different aspects of 
electron-phonon scattering in graphene see Refs. [25l - l4^ . 

A remarkable property of flexural phonons is a cru¬ 
cial role played by anharmonic effects. In particular, 
golden-rule calculation of scattering rate on the deforma¬ 
tion potential created by flexural phonons leads, in the 
harmonic approximation and with electrostatic screening 
taken into account, to values of the electrical conductiv¬ 
ity that are two to three orders of magnitude lower than 
experimentally observed. This drastic discrepancy im¬ 
plies the existence of a certain mechanism that strongly 
suppresses out-of-plane modes. As was demonstrated 
in Ref. taking into account anharmonic interaction 
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between flexural and in-plane phonons dramatically re¬ 
duces the electron-phonon scattering, yielding values of 
the electrical conductivity that are in good agreement 
with experimental findings of Ref. Hence, a compar¬ 
ison of theoretical results with transport measurements 
demonstrates that anharmonic effects are extremely im¬ 
portant in suspended graphene. 

Suppression of scattering because of anharmonicity 
is very favourable for fabrication of ultra-high-mobility 
graphene structures. Further, the anharmonicity gov¬ 
erns lattice thermal transport in suspended graphene, 
which is currently a subject of intense experimental^^— 
and theoretical research (see, e.g. Ref. and references 
therein). Moreover, anharmonicity plays a key role for 
the fundamental issue of the thermodynamic stability 
of graphene membrane. Indeed, due to the softness of 
flexural-phonon modes they might be expected to be very 
efficient in inducing strong thermal out-of-plane fluctu¬ 
ations and thus destroying the membran o^^'^^ by driv¬ 
ing it into the so-called crumpled phase— This question 
was intensively discussed in the literature more than two 
decades ago ^^i^^~ — in connection with biological mem¬ 
branes, polymerized layers and some inorganic surfaces. 
The interest to this topic has been renewed more re- 
cently^^— after discovery of graphene. 

It was found^^— that anharmonic coupling of in-plane 
and out-of-plane phonons stabilizes the membrane for 
sufficiently low temperatures T, so that the membrane 
is in the flat phase at relatively low T and undergoes the 
crumpling transition with increasing T. The main di¬ 
mensionless parameter characterizing the state of the 2D 
membrane is the ratio of the bending rigidity x, to the 
temperature. For graphene this ratio for room tempera¬ 
ture is quite large, xjT ~ 30. This reflects a remarkable 
rigidity of graphene and implies that graphene remains 
in the flat phase up to the temperatures several times 
higher than the room temperature. Moreover, interac¬ 
tion between flexural and in-plane phonons leads to a 
power-law renormalization of the bending rigidity^^i^iS^ 

X ^ x{q) (X q~'^, for g —>■ 0, (3) 

with a certain critical index 77 . Physically, the increase 
of the bending rigidity, Eq. ( 0 , is a manifestation of the 
tendency of the membrane towards the flat phase. 

Development of a controllable analytical treatment of 
renormalization of the bending rigidity and of the crum¬ 
pling transition is not an easy task. The central problem 
is the absence of a small parameter that would control 
the analysis in the physically relevant case of a membrane 
with dimension D = 2 embedded into the space with di¬ 
mension d = 3. For a membrane with arbitrary D and 
d, a systematic treatment turns out to be possible in two 
cases: for 4 — D = e <C 1 and for dc = d — D ^ 1. In 
both limits, there exists a small parameter that controls 
calculations: e in the first case, and 1 /dc in the second 
case. In particular, a theory of crumpling transition for 
D = 2 and dc » 1 was developed in Ref. Issllssl . while a 
renormalization-group (RG) treatment of the membrane 


elastic coefficients in the limit e <C 1 was first discussed 
in Ref. The value of the critical exponent 77 char¬ 

acterising the flat phase was found to be 



for D = 2 and dc > 1 ^*^ and 


12e 

24 dc 


( 4 ) 

( 5 ) 


for e 1 — 

The scaling of the bending rigidity exactly at the crum¬ 
pling transition point is characterised by another crit¬ 
ical exponent rjcr- This exponent determines the frac¬ 
tal (Hausdorff) dimension of the membrane at criticality, 
Dh = 2D/{4 — D — 77cr). It turns out that for D = 2 
and dc ^ 1 the exponent rja coincides with 77 , Eq. ®, 
in the leading order in 1/dc. On the other hand, in the 
vicinity of D = 4 the exponent rjcr was found to be es¬ 
sentially different from 77 , Eq. ®, and scaling as e^, see 
Refs. [sgl^ . 

Eurther approximate calculation schemes (that become 
controllable for dc ^ 1) have been developed, such as the 
self-consistent screening approximation (SCSA)^ and 
the “non-perturbative renormalization group”—. After 
extrapolation to the physical dimensionality, the corre¬ 
sponding results yield 77 = 0.821 and 77 = 0.849, respec¬ 
tively. Clearly, the extrapolation is not controlled para¬ 
metrically; the scattering between the above values may 
serve as a rough estimate of their accuracy. Numerical 
simulations of the problem rave values 77 = 0.60 ± 0.10 
and 77 = 0.72±0.04 (see Ref. [63 and Ref.[6a, respectively). 

In the present paper we explore the interplay of a static 
disorder and anharmonicities in a graphene membrane. 
Some aspects of such a problem have been discussed 
in the literature in a general context of the membrane 
physic a^^i^^'^^ 1 ^^ 1 ^^ . In particular, it was concluded in 
Refs. [6l6l6^ that for D = 4 — e with small e the flat 
phase is stable with respect to various types of disorder. 
On the other hand, it was argued in Ref. l6J that, in the 
leading order in 1/dc, a flat phase of a 2D membrane is 
destroyed already by an infinitely small impurity-induced 
disorder (in-plane quenched random deformations). The 
later result, if applicable in the physical case d = 3, would 
imply an instability of the graphene membrane with re¬ 
spect to an arbitrarily weak disorder. The authors of 
Ref. speculated that high-order terms with respect 
to 1/dc might cure such an instability. In Ref. 1^. influ¬ 
ence of randomness on the flat phase was treated within 
the SCSA, with a focus on a long-range disorder. 

Although effects of disorder in the flat phase (both 
for D = 4 — e and D = 2) were discussed in a num¬ 
ber of publications, the problem of crumpling transi¬ 
tion in a disordered membrane still remains a challenge. 
This problem has a close relation to formation of rip¬ 
ples which are the static out-of-plane random deforma¬ 
tions of membrane. In other words, ripples look like 
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frozen flexural phonons. In graphene, such deformations 
with the height 3 — 10 A and typical spatial scale about 
several hundred angstroms were observed in a number 
of experiments!While several recent theoreti¬ 
cal papers discussed graphene rippleS) ^°i"‘°i'^^i’^^i^^ devel¬ 
opment of a systematic theory of their formation (which 
would explain the mechanism of rippling and predict key 
parameters of the emerging structure) remains a chal¬ 
lenging problem. 

In the present work, we study the crumpling transition 
and the rippling in graphene with a static quenched dis¬ 
order. We perform a detailed RG analysis for the out-of¬ 
plane (“random curvature”) disorder model. Previously, 
the RG treatment of disorder in membrane o^^i^^i^^ has 
been performed only for dimensionality = 4 — e, which 
cannot be directly employed to the physical graphene 
samples. We assume that the membrane dimensional¬ 
ity D is equal (or close) to 2 and use the 1 /dc-expansion 
which allows us to control the theory and to derive cou¬ 
pled RG equations describing a combined flow of and 
effective strength b of the out-of-plane disorder. In this 
way, we establish the phase diagram by determining 
a critical curve b{>c) separating the flat and crumpled 
phases. We demonstrate that, even deep in the flat phase, 
random fluctuation of membrane tension caused by the 
static disorder may strongly change the critical behavior 
of the bending rigidity. We discuss in detail the behavior 
of disorder in the flat phase. 

Far from the critical curve, i.e. deep in the flat 
phase, disorder decays in a power-law way at the scales 
larger than the Ginzburg length which therefore sets a 
characteristic size for the ripples - static frozen out- 
of-plane deformations observed experimentally in sus¬ 
pended graphene. We find that in the limit L —>■ oo 
static and dynamic correlation functions of V/i behave 
as and L~'^, respectively. Hence, ripples are char¬ 
acterized by the anomalous exponent 2?7 in contrast to 
dynamical fluctuations governed by ry. Furthermore, we 
show that, for sufficiently strong disorder, there exists an 
intermediate range of spatial scales where ripples decay 
much slower, with exponent 77 / 4 . One of remarkable re¬ 
sults is that the RG flow of coupling constants may be 
essentially non-monotonous. Specifically, for sufficiently 
strong disorder (close to the critical curve), the bending 
rigidity decreases at the first stage of the renormaliza¬ 
tion, reaches its minimum, and then starts to grow. In 
this near-critical regime the disorder also changes in a 
non-monotonous fashion; it first increases very slowly 
(logarithmically) with the spatial scale L, then reaches 
a maximum at a certain scale and, finally, decreases 
according to a power law at larger scales. In this case, 
the membrane shows fractal properties which imply its 
multiple folding starting form a certain length scale Li 
and flattens at a much larger scale L 2 . 

We also briefly discuss an in-plane disorder and show 
that it is irrelevant in the RG sense, unless its correlation 
function is highly singular at small momenta, i.e., unless 
it is too long-ranged. Thus, the above conclusions remain 


valid also in the presence of in-plane disorder. The in¬ 
plane disorder may, however, affect essentially the bare 
value of the out-of-plane disorder at the atomic scales, 
which serves as a starting point for the RG analysis. 

We conclude the paper by a comparison of our results 
with available experimental data on graphene ripples. 
While the main focus of the paper is on graphene, the de¬ 
veloped theory is quite general and is expected to be ap¬ 
plicable to other crystalline membranes as well. These in¬ 
clude biological membranes like those of red blood cellS)2S 
oxidised graphite or graphene;^ graphane^^ molybde¬ 
num disulphide;^ and boron nitride^S membranes. Fur¬ 
ther examples may likely emerge in near future, in view 
of current active works on engineering of 2D materials 
and structures. 


II. FORMULATION OF THE PROBLEM AND 
MEAN-FIELD ANALYSIS 

To begin with, we note that the energy of membrane 
consists of kinetic and elastic contributions. In this 
paper, we treat the problem semiclassically (see Ap¬ 
pendix 1^ for discussion of the region of applicability of 
the quasiclassical approximation). The kinetic energy 
depends on momenta only. Within the semiclassical ap¬ 
proximation, the phonon momenta can thus be integrated 
out from the very beginning. In what follows we neglect 
the kinetic term (see Refs. for a discussion of 

some effects related to this term) and focus on the study 
of the elastic terms. 

We start from the clean case and consider a general 
Landau-Ginzburg form of the elastic free energy)^ 

F = J ||(a„5„R)2 - ^(a„Ra„R) (6) 

-bM(5aRa^R)2 -b 7;(a„R9cR)^} , 

Here D is the dimension of the membrane, d is the di¬ 
mension of the embedding space and a,/? = 1,..., D. In 
principle, Eq. ([6]) allows one to describe both flat and 
crumpled phases of the membrane. In the latter case, 
one should include an additional term that prevents self¬ 
intersections. In this paper, we focus on the flat phase, so 
that we omit this term. The d—dimensional vector R de¬ 
pends on the point of the D-dimensional reference space, 
X = (cci,..., Xu), i.e., R = R(x) (see Fig. [T]). Equation 
dni) can be obtained by following considerations. Eirst, 
the translational invariance implies that the free energy 
depends on R via the derivatives 9„R only. Second, 
the rotational invariance requires that the energy should 
be scalar with respect to both the embedding (R) space 
and the reference (x) space. Finally, keeping the leading 
terms in the expansion in the field 9a R and in gradients 
yields Eq. dH]). 

Eollowing Ref. [13, we first consider the mean field ap¬ 
proximation which amounts to the linear ansatz, 

R = .^ox. 
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d-dimensional space 


FIG. 1: Membrane with dimension D is embedded into a 
space with dimension d. A point on the membrane surface is 
labeled by a d—dimensional vector R which depends on vector 
X belonging to an arbitrary D—dimensional reference plane. 


This yields the free energy of Landau type, 

F <x-^lt + 2^^{u + Dv). (7) 

Next, we find by minimization of the free energy ©: 

f 77— ^ n X ; for f > 0 

dF/d^o = 0 ^ i (8) 

[ 0, for t < 0 

The case with > 0 corresponds to the flat phase, while 
for ^0=0 fbe manifold R(x) shrinks to a point R = 0, 
which implies that the membrane is crumpled. (Taking 
self-avoidance into account would lead to a hnite radius of 
the crumpled membrane.—) Assuming (in spirit of Lan¬ 
dau theory of phase transitions) that t (x Tc — T, we find 
that oc Tc — T. Hence, the model Eq. m shows a crum¬ 
pling transition at T = already within the mean-field 
approximation. 

We assume now that the system is in the flat phase T < 
Tc and take into account thermal fluctuations around the 
mean-field solution. For this purpose, we write 

R = Cor, (9) 

where 

r = X -I- u + h (10) 

and vectors 

u = (ui, ...,u_d), h= {hi,...,hd-D) 

represents in-plane and out-of-plane displacements, re¬ 
spectively. Substituting Eq. ([S]) into Eq. ([5]) and choosing 
for ^0 the mean-field solution Eq. ([S]), we find 

F = J d^x{|(Ar)2 (11) 

+ ^[dardi3r - - D]'^^ . 

Here 

x = w^l (12) 


is the bending rigidity and 

= 4<;^, A = 8uCo 

are in-plane elastic constants. We see that these con¬ 
stants as well as k turn to zero at the transition point: 
/i, A oc (Tc — T)^, X (xTc — T. 

Next, we use the parametrisation CQD and rewrite 
Eq. (fTTl) in terms of u and h. We assume that spatial 
derivatives of fields are small, |9ah| ^ 1, |(3au| <C 1, so 
that one can neglect terms daudisu compared to i9q,u in 
the second and third contributions to the energy m- 
Further, one may neglect the term >!:(Au)^ compared to 
fi{daU 0 )‘^ and A(i9ctUQ)^ provided that the characteristic 
length scale of the variation of the membrane displace¬ 
ment, q~^, is sufficiently large: 



Under these assumption, one can rewrite Eq. dill) for 
membrane energy in the standard textbook formal: 

F = I |^(Ah)2 -b -b , (14) 

where 

UaP = ^ {daVdpr - 6a0) ^ {OaU/d + dpUa + da^ldph) 

(15) 

is the deformation tensor. 

The following comment is in order here. The question 
about the range of applicability of the simplified model, 
Eq. (imi is in fact somewhat more subtle than one might 
conclude from the above discussion. The point is that the 
elastic constants A, and x are strongly renormalized 
due to anharmonicity and critical fluctuations that it in¬ 
duces. This effect becoming prominent for q <C <?*, where 
q* is the inverse Ginzburg scale, see Eq. (HD) below. At 
such scales, the criterion (IT^ for neglecting the x{Au)‘^ 
term should be modified. In particular, for T = 2, the 
screening by critical fluctuations reduces fj, and A to the 
value ^ U/LIq ^ x^q^jT (see section UlIBI) . Hence, q 
drops out from the estimate Eq. (113 and instead we get 
the condition x: T. As will be demonstrated below, 

the crumpling transition corresponds to x: ~ where 
dc = d — D. Hence, the simplified model (11^1) is suffi¬ 
cient for a description of the crumpling transition pro¬ 
vided that dc ^ 1. On the other hand, it is exactly the 
latter condition that we will use to develop a controllable 
RG treatment for the theory ([HI) . Therefore, neglecting 
the x:(Au)^ term and thus restricting ourselves with the 
free energy (fHl) is fully consistent with the subsequent 
RG analysis. 

Let us now introduce a quenched disorder into the 
model Cl- We will mainly focus on the out-of-plane 
(random curvature) disorder in this paper— Our moti¬ 
vations for analysing the effect of this type of disorder are 
twofold. First, we will show that the random-curvature 
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disorder influences the RG flow (and thus the crumpling 
transition) in an essential way. This should be contrasted 
with the in-plane disorder, which is RG-irrelevant (un¬ 
less its correlation function is highly singular at low mo¬ 
menta), as will be shown below. Second, such disorder 
induces static out-of-plane deformations and is thus di¬ 
rectly related to the effect of rippling. 

The free energy including the random-curvature disor¬ 
der reads^ 

^(u,h) = y {Ah + (3)^ + +. (16) 

Here f3 = /3(x) is a random vector with Gaussian distri¬ 
bution 

^(/3) = exp (^-^J , (17) 

where b is the disorder strength and Zp is a normaliza¬ 
tion factor. In order to study the membrane properties 
in a systematic way, one should take into account fluctu¬ 
ations, thus going beyond the mean-field approximation. 


III. BEYOND MEAN FIELD 

In this Section, we develop an RG treatment of a dis¬ 
ordered membrane. Previously, the curvature disorder 
was treated within the RG approach^ only in dimen¬ 
sionality D = 4 — e. We will see below that the structure 
of the RG near H = 2 is crucially different from the 
case D = A — e. A lowest-order perturbative treatment 
of the curvature disorder in physical dimension D = 2 
was briefly discussed in Ref. roll . However, this consider¬ 
ation is only applicable for short distances (smaller than 
then the Ginzburg length), whereas the most interesting 
physics (captured by the RG) developes on larger scales. 

Our analysis is based on an expansion around the or¬ 
dered (flat-membrane) phase, and thus bears analogy to 
a (T-model treatment of conventional critical phenomena. 
In full similarity with the cr-models, H = 2 is a logarith¬ 
mic dimension for the present problem (which manifests 
in >c/T becoming dimensionless in this case), so that our 
analysis will be appropriate at or near D = 2. This is 
highly favourable, since D = 2 \s exactly the physical 
dimension of a membrane. To keep the theory under 
control, we will also assume a large dimensionality of the 
embedding space, dc 2> 1. 


A. Renormalization of stretching parameter ^ 

Above we rescaled coordinates with the scaling factor 
fo which minimize the free energy within the mean-field 
approximation. Beyond mean field one should take into 
account critical fluctuations. Such fluctuations change 
the optimal scaling factor. To consider this effect, we 
modify Eq. as follows: 

r = ^x-fu-f-h. (18) 


Equation (1181) represents a new rescaling of coordinates. 
The mean-field approximation corresponds to ^ = 1. 
Below, we find that ^ flows away from the mean-field 
value due to the fluctuations. Substituting Eq. (HID into 
Eq. dill) and using the same approximations that were 
done in the course of derivation of Eq. (USD, we arrive at 
the following expression for the free energy: 


F = 


DL^{H + XD/2) 
4 

2(^2 _ 1 ) r 


D 


LD 


- 1 )^ 




(19) 


+ F{u,h), 


where F{u., h) is given by Eq. (IT6l) with u replaced by 

u = ^u. 

Since the product dahdah in Eq. m is averaged over 
the whole volume of the membrane, its fluctuations are 
negligibly small in the limit L —^ oo. Hence, one can 
replace 


dahdah —(clahclQh), 

where (• • •) stands for averaging with the statistical fac¬ 
tor exp[—E(u,h)/T]. Minimizing thus obtained energy 
with respect to we find that the optimal value of ^ 
deviates from the mean field value ^ = 1 due to the fluc¬ 
tuations. 


= 1 - ^ (9„h9„h). (20) 

In order to calculate (clahclah), we notice that E(u, h) 
contains linear, quadratic, cubic, and quartic terms with 
respect to variables h and u. In this section, we neglect 
cubic and quartic terms, thus neglecting anharmonicity. 
The anharmonicity-related effects will be included in the 
next section in the framework of the RG formalism. In 
the harmonic approximation, the u-dependent terms in 
the free energy do not couple with h-dependent ones and, 
therefore, can be integrated out from the very beginning. 
The h-dependent part of the energy in the harmonic ap¬ 
proximation is given by / d^x ^(Ah -)- /3)^/2. An aver¬ 
aging over fluctuation of h and over disorder yields 

aha„h) = d,(I + 6)/^^. (21) 

For D = 2 the integral diverges logarithmically. We will 
see below that U = 2 is the lower critical dimension of 
the theory in the limit of inhnite dimensionality of the 
embedding space d —>• oo. The special role of = 2 is 
not unexpected, since the theory that we are developing 
is based on an expansion near the ordered phase (£ > 
0) and thus bears analogy with cr-models (cf. Ref. 1^ . 
Also, this role of H = 2 can be foreseen already after a 
brief inspection of the free energy (HU), (HZl), since both 
coupling constants T/k and b have a dimensionality of 
L^~'^ and thus become dimensionless in 2D. 

In a vicinity of the logarithmic dimension D = 2, an 
RG formalism can be developed in the conventional way. 
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Introducing an infrared cutoff k, making change of vari- functional integral over u. Hence, we arrive at the free 
ables energy that depends on h vectors onlyi^ 


|2 ^ ^2f.2-D^ 


( 22 ) 


and using Eqs. (EOl) and m, we find an RG equation 
that determines a flow of the stretching parameter ^ with 
the spatial scale k~^, 


dA 


{D - 2)^2 



where A = In (km/k) and 

7 _ 

“ (47r)(^-2)/2r(i:)/2 + l)' 
Exactly at £> = 2, Eq. (|23)) simplifies to 


dA Att \k / 


D = 2. 


(23) 


(24) 


(25) 


The renormalization of is terminated by fc = 1/L. The 
first and second terms in the r.h.s. of Eq. (1251) describe 
contributions of dynamic fluctuations (flexural phonons) 
and of static deformations (ripples), respectively. Equa¬ 
tion (l25|l predicts that ^ becomes zero at a finite system 
size. In other words, a 2D membrane should be crum¬ 
pled in the thermodynamic limit due to both dynamic 
and static deformations. While this conclusion is remi¬ 
niscent of the Mermin-Wagner theorem^ that forbids a 
long-range order in 2D systems, it turns out to be wrong 
(which is a manifestation of the fact that the Mermin- 
Wagner theorem is not applicable to the problem under 
consideration^) . Specifically, inclusion of anharmonicity 
leads to renormalization of k and b and restores the or¬ 
dered (flat) phase in 2D for a finite dimensionality of the 
embedding space, d < oo. 

We are now going to derive RG equations describing 
the renormalization of elastic constants and disorder. We 
will see that at sufficiently large spatial scales, and 
A are screened to values independent on the bare ones. 
The effective theory at such scales is thus determined by 
three running couplings: the stretching parameter (“field 
renormalization”) the bending rigidity (“stiffness”) x 
and the effective disorder strength b. We start from a 
discussion of renormalization of x in the clean case (b = 
0) in Sec. IIII Bl and then derive the full set of RG equation 
for disordered membrane in Sec. IIII Gl . 


B. Renormalization of the bending rigidity in the 
clean case 

We start from the clean-membrane free energy F(u, h). 
Since the theory is quadratic with respect to longitudinal 
modes u, they can be integrated outj ^^'^^ The stretching 
parameter ^ enters E(u,h) via u = only, and evi¬ 
dently drops out after changing variables u — >■ in the 


^ ^ ^ J (26) 

+ ^ J{dkdk'dq)Rq{'k,k') (hk+qh_k) (h_k'-qhk'). 

Here and below we use a short-hand notation (dk) = 
(i^k/(27r)^. The kernel of quartic interaction between 
transverse modes (see Fig. [2|) can be cast in the form 
(see Appendix IbI) 


Rq(k,k') =Ar 


D- 1 


+ M 



D - 1 


(27) 


Here 


= Pk = k - 


q(kq) 


where P is the projection operator related to the trans¬ 
ferred momentum q, 

Pa^ — *^a /3 5 ( 28 ) 


and 


w- _ m(2m + 

{2fi + \)T ’ 



(29) 


Hence, for a generic dimensionality D of the membrane, 
there are two coupling constants, N and Af, controlling 
the interaction strength. An important exception is the 
case D = 2 where the constant M. drops out from the 
theory because the term in the square brackets in Eq. (EZD 
turns to zero. Hence, for D = 2 the interaction reads 


l./^ 2/x(/r -b A) [k X q]^ [k' x q]^ 

(2/i + A) g2 ^ 


(30) 


The bare propagator (which is exact in the absence of 
interaction, A/” = AI = 0) is given by 

{hlht^,) = {27T)^S{k - k') , (31) 


where 




T 

xk'^ 


(32) 


The interaction coupling constants get screened in 
analogy with conventional charges in a media with a fi¬ 
nite polarizability. Evaluating the polarization operator 
to the one-loop order [which is nothing but the random 
phase approximation (RPA)], we find the screened cou¬ 
plings: 


i+Ar(D + i)nq’ 

A. 

1 -b 2Mn^ ■ 


(33) 

(34) 
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k + q -k'-q k + q -k'-q 



FIG. 2: Top: h'^— interaction vertex for generic D (left) and 
for D = 2 (right). Bottom: One-loop diagrams describing 
renormalization of the interaction and of the propagator in a 
clean membrane. 


Thus, the coupling constants M and M become q- 
dependent and are screened independently of each other. 
There is, however, an invariant subspace of the elastic 
coefficients. 


2 fi 






D + 2’ 


(35) 


where the coupling constants stay equal up to a numerical 
coefficient: 




2Mq 
D + 1 ' 


(36) 


The polarization operator Ilq reads 


nq = I(dk)kiGlGl_^. (37) 


to the polarization operator: 


ATq 


1 

(i^ + i)nq’ 


Mq 


1 


(40) 


For the invariant subspace (l35|) . Eqs. (l40l) hold for q <C 
q*, where 




(41) 


is the inverse Ginzburg length which separates the nor¬ 
mal region (q > q*) from the critical one {q < q*). For 
the case when bare couplings M and M are essentially 
different, there are two different scales q^ and q^ at 
which they become screened according to Eq. (HOI) . For 
simplicity, we will assume that the bare values are of the 
same order, M ^ M, so that q^ ^ q^ ^ q*. For a 2D 
membrane this question does not even arise, since the 
coupling A4 is simply absent. 

We are now ready to evaluate the renormalization of 
the bending rigidity by (screened) interaction. The renor¬ 
malized propagator of the /i-field is given by 


^ r 

^ Kk^ + Ek ’ 


(42) 


where the one-loop self-energy reads 


^ 2r ,Arq + (D-2)Mq,4^o,, , 

Ek = — / {dq)^ —- 2.k^G^{k - q). (43) 


For g <C (?*, the effective screened interaction [ACq + {D — 
2)A4q\/(D — 1) that enters Eq. (liSll becomes 


Mq + iD- 2)Mq ^ D 

D-l '^2(D-bl)nq’ 


(44) 


Let us now substitute Eq. (PI) into Eq. (H51) and con¬ 
sider the behavior of Ek at low momenta k q*. Using 
Eq. (1551) for the polarisation operator, we find that the 
integral in Eq. (l43l) scales as fc^ln(l/fc). This implies a 
logarithmic renormalization of the bending rigidity: 


Using (l32|) . we get 


Uq = Ap) 


T 2 


^2q4-D 


where 


(38) 


S>c(k) 


. kjq^-^ 

d,{D + l)AD J ^ '^^fc4|k-q|4 

0 

, (45) 


Y (R)y 

A — \2 J ^ \ 2 J /oq\ 

t,{d-i)/2 22D+1 ^ ’ 

For D = 2, we get A 2 = I/IGtt. 

It follows from Eq. (1551) that for any D < 4 the polar¬ 
ization operator increases according to a power law with 
decreasing q. Hence, as seen from Eqs. dSSl) and dMl, 
for sufficiently small q couplings M and M. become inde¬ 
pendent of their bare values and inversely proportional 


where we neglected the contribution from q > q* which is 
convergent for fc —> 0 and is therefore small. The constant 
r] can be expressed as 


7 ] = 


1 D r H(£),7?o)?7o 
dcD + l _ A{D, 770 ) 




(46) 


with integrals A{D,ri) and B{D,r]) as defined in Ap- 
pendix[C] The functions A{D, rj) and B{D, rj) will be used 



























below for comparison with SCSA<^ Note that A{D, 0) = 
Ad- 

For a sufficiently large spatial scale k~^ the correction 
(H51) ceases to be small and gets promoted to a RG equa¬ 
tion, 


d>c 

'dK 


= 


(47) 


where A = \n{q*/k). We thus see that rj is the anomalous 
dimension of the bending rigidity. Using Eqs. (ED and 
(IC2I) . we find 


, mD-D2‘> _ Tj^ _ 

r(f)r(2-f)r(2 + f)- 

It is worth pointing out that, at variance with Eq. (ED, 
the integral (I45|) is always logarithmic (i.e., also away 
from D = 2). Thus, Eq. (H51) is applicable for any D 
as well. The only condition of validity of this equation 
is dc ^ 1) which implies that 77 <C 1. The equation 
simplifies for I? = 2 and D = 4 — e, yielding 

f 2/4, for D = 2, 

(49) 

1 12e/dc, for D = 4 — e. 


[Since dc 2 > 1, the second line in this equation agrees with 
Eq. (ID) obtained by e—expansion.] For generic values of d 
and Z?, the exponent 77 is not small. To find it, one would 
have to take into account higher-loop contributions to 
RG equations. Such a calculation does not appear feasi¬ 
ble because of the absence of a small parameter. On the 
other hand, one may develop a self-consistent extension 
of the one-loop theory by inserting renormalized Green 
functions into the one-loop diagrams, which amounts to 
the replacement^^ of with G in Eqs. (1571) and (H51) . 
Corresponding calculations yield a self-consistency equa¬ 
tion for 77 


1 D 

d^D + l A{D, 77 ) ’ 


(50) 


which was derived in Ref. 1^ . This equation can be ob¬ 
tained from Eq. (H51) by replacing 770 with 77 . Although 
such a procedure is not controlled parametrically for 
physical membranes (d = 3, Z? = 2), it gives in this case 
a value 77 « 0.821 which turns out to be in a reasonable 
agreement^! with numerical values 77 = 0.60 ± 0.10 and 
77 = 0.72 ± 0.04 (see Ref. HI and Ref. HI, respectively). 


C. Renormalization group for disordered 
membrane 


Now we include in the consideration the random cur¬ 
vature disorder which modifies only the first term in 
Eq. (1261) . In the coordinate representation, this term be¬ 
comes 



(51) 


8 

where /3 = /3(x) is a random vector with Gaussian distri¬ 
bution ED. Disorder averaging can be performed with 
making use of the replica trick. To this end, we intro¬ 
duce N replicas of the field h (i.e., make a replacement 
h —>■ h" with a replica index n= 1,... ,N) and replicate 
the free energy F : 


^rep ^ ^rep ^rep^ 


(52) 


where 


n^N 


9 {dk)k%Z+f3^f 


(53) 


and 


n^N 


.o — .y - 


n—1 


(54) 


X / {dkdk'dq)R^{k,k’) (hD+qh((k) (h!!k'-qK,) . 


Next, we average exp(—F^ep/T) with P{I3), thus arriving 
at the following effective action: 

^ E / + Ff P, (55) 

n,m 


where we have introduced a replica-space matrix k with 
elements given by 


yi = K — 


byp' 

T + bxN 


J. 


(56) 


Here J is a matrix with all elements equal to unity: 
jnm _ bare propagator thus becomes a matrix 

in the replica space: 

where 

/=^ (58) 

is the dimensionless parameter characterising a ratio be¬ 
tween the bare disorder, b, and the bare dynamical fluc¬ 
tuations, T/k. 

The polarization operator also becomes a replica-space 
matrix. In the one-loop order (i.e., within RPA), its ele¬ 
ments read 


yrnm _ 

q “ £)2 _ 


^ /(dfc)fclGr™G°'7r 


(59) 


F2 - 1 


{dk)F 


(A-i)"™ (A-i)"™ 

■ P |q-k|4 ■ 


Using the property (5"™ -b = (1 + 2.f)5^^ + 

/V"™, we hnd 


Hq = Ad 


T2 




1 + 2/ + /V 


(Ah -f (3)^ , 


(60) 
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Substituting Eq. (l60l) into Eqs. (j33ll nnd (13411 . wg obtuin 
the screened couplings AAq and as matrices in the 
replica space. The Ginzburg scale q* is now affected by 
the disorder that enters the polarization operator (1501) : 


mT(1 + 2/) 


1 i/(4-r') 


(61) 


It is worth noting that for strong disorder or low tem¬ 
peratures (/ 1) g* ~ is independent of 

temperature, while for weak disorder (/ -C 1) we recover 
Eq. dH]), q* oc 

Let us now calculate the self-energy determining the 
renormalization of the bending rigidity. In similarity with 
Eq. (l4^ , the effective interaction U with matrix elements 
Unm (see Fig. [3]) is expressed for g <C g* in terms of the 
polarization operator, 



U = 


2(D-f 1) 

D l + 2f + pN - pj 

2{D + \) T^Ad (l + 2/)(l + 2/ + /2iV) 


(62) 


The replica generalization of Eq. (l43l) for self-energy thus 
takes the form 


S nm 

k 


2T 

dc 



^ ) ^0,nm 

2(D + 1) *^-'1 


(63) 


Substituting here the Green function, Eq. dsm, using the 
property 


1 + 2/+ pN-pj (i + fjy 


l + 3f + p(N + l) + pN-pJ 


and taking the limit IV 0, we finally arrive at the 
matrix equation which governs the renormalization of >c 
in the presence of disorder: 


dx 1+3/ + /2-/3J 

dA (l + 2/)2 


(64) 


where rj is given by Eq. (|46|) . Substituting Eq. (l56|) into 
l.h.s. of Eq. (l64l) . and separating terms proportional to J 
from scalar ones, we find two equations describing renor¬ 
malization of the bending rigidity and of the disorder 
strength: 


dx 

1 + 3 / + P 

(65) 

dA 

(l + 2/)2 ’ 

d{x/) 

P 

(66) 

dA 

'^''{1 + 2 /y 


Equations (l23l) . (l65l) . (l66l) represent the full set of RG 
equations which determine the behavior of the system 
for g ^ g* [the Ginzburg scale g* is defined here with 
the bare disorder strength /o in Eq. mi To analyse 


FIG. 3: Replica-index structure of the interaction and 
of the one-loop diagrams for a disordered membrane. 


the fixed points and the scaling flow, we rewrite them 
in terms of appropriate dimensionless couplings. To this 
end, we introduce a rescaled bending rigidity, 

k = Kp = . (67) 


The meaning of Eq. [53 is twofold. First, the ratio k/T 
is made dimensionless (also for a dimensionality devi¬ 
ating from D = 2). Second, it takes into account the 
“compression” of the membrane controlled by the fac¬ 
tor ^ in Eq. m, which can be viewed as an analog of 
the order-parameter field renormalization in the conven¬ 
tional cr-model RG. Expressing further the flow of / from 
Eqs. ((551) and (1551) . we cast the RG equation into the fol¬ 
lowing form: 


dA 

dk 

dA 

d/_ 

dA 


PT 

47r k 


(l + /)e^ 


X 


D — 2 T] 


1 + 3/ + /^ ' 

(I + 2/)2 _ 


/(I+ 3/) 
^I + 2/)2' 


( 68 ) 

^(l + /)(69) 


(70) 


It is worth discussing a structure of the derived RG 
equations. Equations (1591) and (1751) describe evolution of 
two dimensionless couplings of the theory: the inverse di¬ 
mensionless bending rigidity T/k and the dimensionless 
disorder b = /T/k. The solution of one-loop equations 
is simplified by the fact that the equation for / (ratio of 
disorder to the inverse bending rigidity), Eq. (iTOl) . decou¬ 
ples and thus can be straightforwardly integrated. After 
this, Eq. (l69l) for k can be solved. Finally, Eq. (1581) is 
analogous to the wave-function renormalization equation 
of the cr-model RG. It is a “slave” equation which can be 
solved at the last stage. 

We emphasize that the set of the RG equations 
dMl), dMl), and (1751) strongly differs from the RG equa- 
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previously derived in dimensionality D = 
4 — e. Apart from the dimensionality (that allows us 
to employ the RG approach to physical D = 2 graphene 
samples), the crucial difference is in the role of screening 
near two dimensions. In our case, the screening of anhar- 
monic couplings eliminates the separate renormalization 
of elastic constants A and /i; they now only determine 
the starting point of the RG, q*, Eq. (EH). For the same 
reason, only the out-of-plane (curvature) short-range dis¬ 
order survives the renormalization, see discussion of the 
in-plane disorder in Sec. IIIIE)] below. 

One can further simplify Eqs. (1551) and (1551) by intro¬ 
ducing a dimensionless coupling 


47r?7i< 
dcTil + f) 


(71) 


and a is the effective disorder strength. Replicating the 
field h and averaging over disorder, we find that the bare 
interaction 17°^ (see Fig. [3|) acquires now off-diagonal 
elements in the replica space: 

J. (78) 

Next, we find the screened interaction by solving the 
equation U = Uq — UqUU (see Fig. |5]), which yields 

U={1 + UoUy^Uo = (1 + (79) 

For q q*, the polarization operator grows as 1/g^, so 

that one can neglect the term 11“^ As a result, 

both the elastic constant Af and the disorder strength a 
drop out from the effective interaction: 


[for D = 2 this transformation takes the form K = 
(87r/dgT)>ir/(l-|-/)]. We also change the variable A to 


z = ryA, 

thus arriving at the following RG equations: 

de e 


dz 

dK 

dz 

dz 


77’ 


= K 


£2 


1 + 5/ + 7/2 + f 
(l + /)(l + 2/)2 

/(I+ 3/) 

(l + 2/)2- 


Here 


£2 = 


D-2 {D- 2)da 




- 1 , 


(72) 

(73) 

(74) 

(75) 

(76) 


Before proceeding to a detailed analysis of these RG 
equations, we briefly discuss the effect of in-plane disor¬ 
der. 


D. In-plane disorder 


In most of this paper, we explore the problem in the 
presence of out-of-plane (“random curvature”) disorder. 
In this subsection we briefly analyse the in-plane disorder 
and show that, with an exception of the case of topologi¬ 
cal defects (disclinations), it can be safely neglected. We 
restrict ourselves to the case D = 2. Let us consider an 
impurity-induced isotropic in-plane disorder that leads^ 
to the following modification of the free energy (1551) : 


= (dfc)fc^lhkP 


(77) 


(k X q)2 


Ik+q 


h-k + Cc 


Here Cq is the Fourier transform of the random quenched 
field c(x) with the Gaussian distribution 


P[c(x)] oc exp 



U^U-\ tor q<^q*. (80) 

We thus conclude that, in the critical region g <C g*, the 
in-plane random potential becomes irrelevant and does 
not affect critical behavior of >c. 

This conclusion should be contrasted with that of 
Ref. 113 where it was argued that an arbitrarily weak in¬ 
plane impurity disorder destroys the flat phase. Specif¬ 
ically, it was found that the membrane is destroyed for 
a > (Tc*^^, where 1 /is strongly divergent in the in¬ 
frared limit: cx Here, gm is the infrared 

cutoff momentum (inverse system size, giR ^ 1/A). One 
can demonstrate that in fact screening cures this diver¬ 
gency (cf. Ref. 113 ) and leads to a replacement of giR 
with g*, thus yielding a finite value CTc oc T/^. A not too 
strong in-plane disorder, (Tq < cTc will thus simply lead 
to ultraviolet corrections to the bare parameters (out-of- 
plane disorder bo and bending rigidity xtq) of our the¬ 
ory. Specifically, a natural expectation is that the in¬ 
plane disorder would enhance bo and reduce xq . In other 
words, the in-plane disorder is RG-irrelevant at length 
scales > 1/g* but its contribution on shorter (atomic) 
scales can affect the starting point of RG. 

The above result on irrelevance of the in-plane disor¬ 
der was based on an assumption of a finite-range disorder 
c(x), which means that the disorder strength Uq remains 
finite in the limit q —>■ 0. It is easy to see that this 
assumption can be in fact weakened: it is sufficient to 
require that the singularity of Cq at q —>■ 0 (if present) 
is not stronger than l/g2““ with a > 0. This corre¬ 
sponds to spatial correlations with a power-law decay, 
(c(O)c(x)) oc x~°‘. In this case, inverting Eq. (175)) and 
taking the limit TV —>■ 0, we find 


U-\N^0) = ^ + aJ. (81) 

Since —)► g“ for g —>• 0 , we conclude that ?7g"^n“^ 

vanishes at small g, so that the derivation of (I50l) retains 
its validity. Only for a —> 0 (which corresponds to log¬ 
arithmic real-space correlations) we get crqll~^ —>■ const 
for g —>■ 0. Such a disorder is RG-marginal and thus 
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will influence the flow of other couplings. Physically, the 
1 /q^ in-plane disorder corresponds to random topologi¬ 
cal defects—disclinations (see Ref. for a discussion of 
different types of long-range disorder and their treatment 
within the SCSA). A detailed analysis of this type of in¬ 
plane disorder is outside of the scope of this paper and 
will be presented elsewhere. 


IV. ANALYSIS OF RG EQUATIONS: 
CRUMPLING TRANSITION 


Hence, the fixed point (l86l) separates the crumpled and 
flat phases 

For a membrane dimensionality D slightly deviating 
from 2 , we get instead of (IMll 

(UK 

— =K(e 2 + l)-l. (90) 

az 

Equation (|90l) implies that the lower critical dimension 
for crumpling transition can be found from the condition 
62 = — 1 , which yields 


In this section we analyze the RG equations derived 
above for clean and disordered cases. 


A. Clean membrane 


In the absence of disorder (/ = 0) and for £> = 2, 
Eqs. d73l) and (1741) take the form 


dz 

dK 

dz 


Jl 

K’ 

K-1 


(82) 

(83) 


Here we have taken into account that ^ ^ for H = 2. 

The initial condition for Eq. (15^ is2^ 


^0 = 1, for z = 0. (84) 

From Eqs. (l82l) and (l83l) we conclude that there exists an 
unstable fixed point 


iFcr = l 


(85) 


or, equivalently, 


U„ = 2-^, (91) 

dc 

in agreement with previous studiesi^. 

Exactly at the transition point, when K = K^r = 
1/(1- 1 - 62 ), the stretching factor ^ decays with L according 
to a power law: 


(92) 

where t = {D — 2 + rj)/2. In other words, the extension 
of the membrane R in the embedding space scales with 
its “intrinsic” length L as R = The exponent r 

determines thus the fractal (Hausdorff) dimension Djj 
of the membrane at criticality (defined by the relation 
L^), yielding Dh = D/{1 — t) = 2D/{A —D— rj). 
Let us discuss in more detail geometric properties of 
the membrane, which are determined by the behavior of 
To this end, let us consider two points ri and r 2 on 
the membrane and the corresponding points xi and X 2 
in the reference plane. According to Eq. (ITSl) and to the 
RG procedure, we have a scaling relation 

(|ri - r 2 |^) ~ (xi - X 2 )^ (93) 


^cr = (86) 

Indeed, assuming that the starting value of the bending 
rigidity, xq = >Cz=o = Xz=o, exceeds the critical value, 
xq > Xcr, we find from Eqs. (1821) and (l83l) 




>^cre ^ + Xq — X^ 
Xo 


and, consequently, 


^2 -^0 -^ci 

S2=00 ~ 

Xo 


(87) 


( 88 ) 


Therefore, above the critical point membrane remains in 
the flat phase in the course of renormalization. On the 
other hand, one can easily check that below the criti¬ 
cal point (for xq < x^r), the membrane shrinks to the 
crumpled phase, ^ = 0 , at a finite scale 


L 


q 


-1 


1 / ^CT 
q* \ Xcr - Xo 


de/2 


(89) 


where is the value of the running renormalization- 
parameter ^ at the RG spatial scale x. When ^x drops 
down to a value substantially smaller than unity (say, to 
1/2), the embedding-space distance |ri — r 2 | between the 
points become essentially smaller than the intrinsic dis¬ 
tance |xi — X 2 I. This indicates that a membrane of such 
size starts to show “folding” (i.e., strong spatial variation 
of the normal vector to the membrane in the embedding 
space). 

Deep in the flat phase, the renormalization of ^ at the 
whole interval of RG scales, l/q* < L < 00 , is relatively 
weak. This implies that the membrane does not fold. In 
other words, although the surface of a membrane is not 
exactly flat due to dynamical fluctuations (and also due 
to ripples in the disordered case, as discussed below), the 
spatial variation of the normal vector remains relatively 
small. 

On the other hand, when xq approaches the critical 
value, Eq. (l86l) (i.e., when xq — Xci- >Ccr), the mem¬ 
brane shows fractal folding at the broad interval of length 
scales. Specifically, as seen from Eq, (1571) . this folding be¬ 
comes strong for z ^ 1, i.e. at the length scale L ^ Li, 
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where 


Li ^ 



At a much larger scale, 


L 2 ^ 


1 

q* 


Kq - Kc- 


1/r] 


(94) 


(95) 


exp(—z) becomes on the order of {>cq — >Cc-c)l>Cc-[ and ^ 
saturates. Thus, the membrane has a fractal geometry 
in the interval Li <C A <C ^ 2 - At larger scales, L ^ L 2 
the membrane flattens. Exactly at the transition, L 2 
diverges and the membrane remains fractal at arbitrarily 
large scales. 


B. Disordered membrane 

1. RG flow in disordered case 


Let us now consider the disordered case. From now on 
and till the end of the paper, we will assume that D = 2 
and, consequently, €2 = 0. We see from Eq. dZSD that 
/ monotonously decreases in course of renormalization. 
Dividing Eq. dm) by Eq. (I75|), we find an equation that 
determines variation of K with /. Its solution yields the 
function K (/) which can be written in the following form: 

Kif) = K^rif) + ^ -< 5 . ( 96 ) 

Here K^Af) is the critical curve in the (AT, /) plane which 
separates the crumpled and flat phases (see Fig.|l|): 


K,Af) 


(3/ + l)i/9e-//3 I ^l + 2y)^(l+y)ey/^ 
/(! + /) (3y +1)10/9 

0 


(97) 

The phase boundary has the following asymptotic be¬ 
haviour at small and large /: 


Kr 


1 + /^/3, for / -t 0; 
4(1-5//), for/-)- 00 . 


For / = 0 we recover Eq. ((551) as expected. The param¬ 
eter S controlling the deviation of the RG flow line from 
the critical one depends on initial values of / and K in 
the following way: 

The flat phase corresponds to <5 > 0, while the crumpled 
phase to ^ < 0. 

The RG flow in the (AT, /) plane is illustrated in 
Fig. dH). One observes that in the flat phase the evolution 
of AT is non-monotonous. For ^ <C 1 the minimum of AT 


K 




FIG. 4: Upper panel: Critical curve (<5 = 0) separating the 
crumpled and flat phases in the K - f plane. Lower panel: 
RG flow lines in the K - f plane for a disordered membrane 
for different values of S (S= -200, -100, -50, -3, -0.4, 0, 1.2, 
7.5, 50, 100, 200 from bottom to top). 


lies in the region of small /, where AT Ri l-|-/9/3-|-i5//. 
Minimization yields the position of the minimum: 

Un - (3<5/2)i/3, « 1 + (3(5/2)9/3, for 5 < 1. 

( 100 ) 

In the opposite case, 5 1, the minimum corresponds 

to large /, where K Ri 4—20//-|-95exp(—//3)/(3/)^^/9. 
We find with logarithmical precision 

Uin-31n(S), K^in-4-20/31n(S), for ,5 » 1. (101) 


2. Disorder-induced crumpling 


Let us fix initial value of >!ro = = >*^ 2=0 and con¬ 

sider what happens with increasing disorder strength b. 
Since /o = Ij we have the following initial values of AT 
and /: 


Ko = 


^0 

>rcr(l + /o) 


/o 


bo>io 

T 


( 102 ) 


where xTcr is given by Eq. (1861) . The crumpling occurs 
when Kq becomes smaller than A'cr(/o)- Hence, the criti¬ 
cal curve in the (xtoi ^o) plane is implicitly determined by 
equations >ro = >fcr(l + /o)A'cr(/o) and bo = /oT/xq = 
{87r/(f)fo/ (l-|-/o)A:cr(/o), with /o varying in the interval 
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bpdl 



FIG. 5: Crumpling-transition phase diagram in the plane of 
initial parameters (bending rigidity >co and disorder 6o). 


(0,oo). The dependence bo{^o} found from these equa¬ 
tions is plotted in Fig. ©■ The asymptotic behaviour of 
the critical line reads 

, _ ^ f 4(>ro - >^cr)/>fcr, for xro >fcr; 

^ ~ \ 1 + 19>rcr/16>^o, for xq ^ oo. ^ ' 

We notice that the dependence is non-monotonous. In 
other words, there exists an optimal value of bare bending 
rigidity, xtq ~ 4.2x:cr, for which the membrane is most 
robust to disorder-induced crumpling. It is also worth 
stressing that the critical disorder saturates in the limit 
of infinite bending rigidity at the universal value bp = 
27i-/d2. 

V. RIPPLING AND FOLDING OF MEMBRANE 
IN THE FLAT PHASE 

We have thus established the phase diagram of the 
crumpling transition in the plane of initial parameters, 
i.e., bending rigidity xtq and disorder bp. An important 
question that remains to be explored is the evolution of 
physical observables with the length scale. This will be 
done in the present section. We will, in particular, show 
that the disorder strength decreases beyond a certain 
length scale L^. 

After rescaling the disorder strength 


bd^ 



FIG. 6: RG flow in the (it,6)—plane for S = -200, -100, - 
50, -3, -0.4, -0.3, -0.2, -0.1, 0, 0.1, 0.3, 1.2, 7.5, 50, 100, 200, 
increasing from top to bottom. 


Fig. [S] The RG flow lines are illustrated in Fig. [51 We 
see that in the flat phase at the first stage of the renor¬ 
malization the bending rigidity strongly decreases, while 
the disorder strength slightly increases. This behavior 
indicates the tendency of the membrane to the disorder- 
induced crumpling. However, below the critical line the 
disorder strength is insufficient to destroy the membrane, 
so that xr eventually shows a minimum and then starts 
to grow, while the disorder gets suppressed. Close to the 
minimum of xr, the RG curve is approximately vertical. 
This behavior can be interpreted as screening of disorder 
by softened membrane. 

The minimal value of the bending rigidity is given by 


r 1 -b 2-\A, for (5 —>■ 0, 
\ 12 In 5, for (5 —>■ oo. 


(107) 


A particularly interesting behaviour is predicted slightly 
above the critical curve (<5 < 0 with |(5| <C 1). In this 
case, the disorder first increases, then reaches a max¬ 
imum, starts to decrease and, finally, after reaching a 
deep minimum, increases again and goes to infinity. 

The geometry of membrane is determined by the be¬ 
havior of the stretching factor which can be expressed 
in terms of K{f) by using Eqs. (1731) and d75|) : 


b = ^, (104) 

one finds that the RG equations in the (iir, &)—plane are 
implicitly determined by the following equations 



— = {l + f)K{f), 

(105) 


bdl 4/ 

(106) 


2^ {l + f)K{fy 

where K{f) is 
interval [0, oo]. 

given by Eq. (ITOl) and / varies in the 
For (5 = 0 we reproduce the critical 


curve bp (tip) separating the crumpled and flat phases, see 


= exp 



(l + 2pd 
y(l + 3y)K(y) 


(108) 


The value of ^ in the limit L —>■ oo is found from this 
equation by putting / = 0 : ^z=oo = C/=o- 

Having in mind applications to graphene, let us now 
discuss in more detail the scale dependence of the bend¬ 
ing rigidity and disorder in the flat phase (5 > 0). Similar 
to the clean case, the behavior of membrane is essentially 
different deep in the flat phase (when S is large) and in 
the near-critical regime (when 6 is small). The regions 
corresponding to different regimes are shown schemati¬ 
cally in Fig. Hon the {S,fp) plane. Regions (I), (H), and 
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FIG. 7: Phase diagram in the {S, fo) plane. Critical curve 
(shown in red) separates flat and crumpled phases. Clean case 
corresponds to horizontal axis {fo ~ 0). Regions (I), (II), (III) 
correspond to a near-critical regime within the flat phase. In 
this part of the phase space, the membrane shows critical 
(fractal) folding at intermediate scales before flattening at 
larger scales. Regions (IV) and (V) correspond to a rippled 
membrane deep in the flat phase. Blue curves correspond to 
fixed values of bare bending rigidity i<o- (a) {^o — Xcr)/>fcr 
1; (b)(x:o — Xcr)/>rcr S> 1. Bare disorder increases along these 
curves from the bottom to the top. 


(II) are sufficiently close to critical line, so that the mem¬ 
brane undergoes strong folding in course of the renormal¬ 
ization. Therefore, we term these regions near-critical. 
In contrast, for membranes with a starting point within 
regions (IV) and (V), the stretching parameter ^ does 
not change essentially in course of renormalization (i.e., 
^ remains close to unity in the limit L —>■ oo). Such 
membranes do not fold and only show small dynamical 
wrinkling and static rippling. 

We consider separately the cases of small and large 
deviations from the critical transition line. 


A. Close to critical line [regions (I), (II), and (III) 
in Fig. (3- 

When the system is in the flat phase but not far from 
the transition (i.e, <) > 0 but i5 <C 1, regions (I) and (II) in 
Fig- El), the dependence i<(/) has a minimum at / ~ \/S. 
The minimal value is given by (fi/Xcr)min — 1 + 2-\/5 and 
therefore is very close to the critical value for the clean 
membrane. The disorder strength shows a maximum at 
a much larger values of /, namely, / ~ 2. The maximal 
value of disorder is given by ~ 1.9. 

For strong disorder, when the starting values of / is 
large, /o ^ 1 [region (II)], the overall RG evolution is 
particularly rich. At the initial stage of renormalization 
(i.e., as long as / remains large), we find form Eq. (TTSl) 
that / changes linearly with 2 : 

fiz) = fo-^, for/»l. (109) 

As follows from Eq. (I106L at such spatial scales the dis¬ 


order slowly increases, 

27r ^ f{z) 


k 


( 110 ) 


Here we have taken into account that (5 <C 1, replaced 
K with Kcr in the denominator of Eq. ()106L and used 
the large-/ asymptotic of the function Ka [see Eq. d^ ]. 
The bending rigidity decreases linearly with z: 


>c 

Xcr 


4/(^) 


( 111 ) 


The scaling of disorder changes when / becomes of the 
order of unity. This happens for 

z,«^ + 0[ln/o] (112) 

As seen from Eq. cni), for large scales / decays exponen¬ 
tially with z, i.e., /(z) ~ exp(—z), so that the disorder 
starts to decrease as 


b{z)dl^^ 4/(z) 
27r ^ K[f{z)y 


(113) 


Although the disorder strength starts to fall, the bending 
rigidity is still strongly affected by disorder and continues 
to decrease up to a value >tmin where / becomes quite 
small (of the order of V^). Only after this the bending 
rigidity begins to grow. 

In the case when the starting value of disorder is 
weaker, <C /o “C 1 [region (I)], it does not show a 
maximum and monotonously decreases in course of the 
renormalization. However, as was mentioned in the pre¬ 
vious paragraph, the disorder strongly affects as long 
as / > Vs. 

The disorder becomes irrelevant for the evolution of 
bending rigidity for fo -C VS. In this case, / ~ 
/oexp(—z) and, using Eqs. (IMI) . (IMI) . (l99l) . (11051) . and 
(I106L we arrive at the following scaling dependencies 


k Ki >fcr + {Xo - >Ccr)ey 

~ boKoe~^ 

b Ri - - --—. 

>fcr + (>f0 - >Ccr)e^ 


(114) 

(115) 


The disorder strength decreases in course of the renor¬ 
malization, first as exp(—z) and, at larger scales, when 
bending rigidity becomes large enough, as exp(—2z). A 
crossover between these two regimes happens at / ~ (5, 
i.e. on the boundary separating the near-critical and fiat 
phases. 

Using Eq. (11081) . one can easily find that in regions (I) 
and (H) the membrane folds in the course of the renor¬ 
malization. In particular, in the region (H), we find that 
at the early stage of renormalization (when / ^ 1) the 
rescaling parameter ^ is given by V ~ exp[(/o — /)/3] ~ 
exp[—z/4] = {q*L)~^/^. Hence, first folding of membrane 
occurs at z Ri 4. This yields the following estimate for the 
length scale of folding: 



(116) 


















15 


One can also find that ^ saturates at an exponentially 
small value 

(117) 

for L ~ L 2 , where 

^2-^ ^ (^y) (118) 

It is worth stressing that L 2 coincides (up to a prefac¬ 
tor) with a length scale corresponding to maximum of 
disorder. 

In the region (I) disorder does not essentially affect the 
folding of the membrane, i.e., the dynamical fluctuations 
dominate. Therefore, Li and L 2 are given by Eqs. ([M|) 
and (|95|l derived above for the clean case. 

Let us now discuss the region (III), where /o 
31ni5 I. Using the asymptotic of K{f) at large / [see 
Eqs. (IMl) and (1551) ]. we find that the functions >c{f) 
and b{f) have, respectively, minimum and maximum at 
fm ~ 3 In d. The minimal value of rigidity is much larger 
than the critical value, i<min/-^cr ~ 12 In 6, while the max¬ 
imal value of disorder slightly exceeds the universal value 
27r/d^; specifically, bd‘^j2Tr — 1 ~ 1/fm- 

In view of the assumption /o > fm, the maximum in b 
is reached at 

= ^(/o - Sind). (II9) 

Eor z Zr the behavior of b and can be well described 
by Eqs. (IllOl) and (11111) with f{z) given by Eq. (I109D . 
Eor z ^ Zr, one can neglect Ka in expression (l96)) . so 
that Eqs. (I105D and (I106D become 


>Ccr 

/ 

(120) 

bdl 

4/^ ,f/3 

(121) 

27T 

<5(3/+ 1)1/9 ■ 


Since /m ^ 1, there exists a large interval of z, where 
these formulas can be used and at the same time / ^ 1, 
so that Eq. (|109p is applicable. In this interval we get, 
with exponential precision, 

k oc b oc (122) 

Finally, for very large z, when / becomes smaller than 
unity and decays as exp(—z), we find from Eqs. (11201) 
and (|121D : 

iiroce^, b oc . (123) 

The latter regime is well described by Eqs. (I114p and 
(IllSP with XTq ~ ^ X^cr- 

In analogy with the regions (I) and (II), we find that in 
the region (III) the membrane undergoes a critical folding 
in a broad range of length scales, Li < L < L2. The value 


of Li is given by the same equation (11161) as in the region 
(II). The stretching factor saturates at the value 

= de-*/3 (124) 

when the spatial scale L becomes of the order of 

L 2 - —e^(*-31n5)/3,)^ (125) 

q* ^ ' 

Again, up to a prefactor this length coincides with the 
length Lt where disorder has a maximum. 


B. Away from the critical line [regions (IV) and 
(V) in Fig. 0. 

The evolution of bending rigidity and disorder for a 
membrane with initial couplings within the region (IV) 
or (V) coincides with the intermediate (or, respectively, 
final) stage of renormalization of k and b for the case 
of bare couplings in the region (III). Specifically, both k 
and b change monotonously with increasing spatial scale: 
k increases, while b decreases. For the region (IV) there 
exists a large interval of length scales L with 1 < / < 
3In(5, where Eqs. (11201) . (11211) . and (11221) apply. With 
further increase of L, we enter the region / < 1, where 
scaling dependencies change to Eq. (I123D . For the region 
(V) the asymptotic formulas (11231) are valid from the very 
beginning of the renormalization. 

Let us analyze the evolution of the stretching factor 
^ in these two regions. For region (V) the membrane 
stretching is fully controlled by dynamical fluctuations, 
while the disorder is irrelevant. Hence, one can use 
Eqs. (1571) and (1551) derived for the clean case. One can 
easily check that within the region (V) xq —>ircr ^ xTcr, so 
that ^ does not change essentially and, consequently, the 
membrane does not fold, remaining approximately flat 
at all scales. The same statement is valid for region (IV) 
as well. Indeed, in this case, one can use the inequality 
1 < / < 31nd and replace the function K{f) in Eq. p08l) 
with /-i^/®exp(-//3) > 1 [see Eq. ^]. Esti¬ 

mating then the integral in Eq. (llOSp . we find that the 
renormalization of ^ remains small. 

This completes the analysis of the RG flow of the bend¬ 
ing rigidity and the disorder for membranes with bare 
couplings in different regions of the flat phase. We have 
seen that the evolution can be very non-trivial, with sev¬ 
eral intermediate scaling regimes. It is worth emphasis¬ 
ing, however, that the scaling behaviour at longest scales, 
L ^ q~^ —>■ 00 , is the same for the whole flat phase: the 
bending rigidity increases, while the disorder gets sup¬ 
pressed according to universal power laws, 

k(xL'^, b (X for L—>- 00 , (126) 

with 7] ~ 2/dc in the large-dc limit. 
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FIG. 8: Scale dependence of the renormalized bending rigid¬ 
ity iir in the fiat phase for <5=1, 5, 25, 100, 300, 1000, 10000, 
100000, increasing from bottom to top. Critical curve {S — 0) 
is shown in red. 


C. Spatial scale of ripples 

As follows from the above discussion, the scaling de¬ 
pendence of both k and b in the flat phase are especially 
non-trivial when the bare disorder is sufficiently strong. 
In particular, b shows in this case a maximum in the 
course of renormalization. The spatial scale correspond¬ 
ing to the maximum can be easily found from expression 
for Zr and is given by 

^ ^ for d <C 1, /o > 1; 

I gd„(/o-31n5)/6^ for(5»l, /o>lnA 

(127) 

For L < Lr the disorder slowly increases with z, while 
for L > Lr it decays exponentially with 2 (i.e., according 
to a power law with respect to L). In other words, the 
disorder is, roughly speaking, “switched off” at L > Lr- 
Hence, one can interpret Lr as a characteristic scale of 
random static deformations—ripples. It is worth empha¬ 
sising at this point that in the case of a nearly critical 
membrane such “ripples” are multiply folded (fractal) 
configurations. 

When the bare disorder is weaker [regions (IV) and (V) 
in Fig. [7], the ripples take a more conventional form of 
relatively small static deformations of a nearly flat sur¬ 
face. Indeed, in this case ^ « I, so that the membrane 
does not fold. For both regions, the suppression of dis¬ 
order begins already on the Ginzburg length 1/q*. The 
scaling of disorder at the initial stage of renormalization 
is different for regions (IV) and (V): 

6 oc (g*T)■’'/^ (128) 

for region (IV) and 

6oc (g*T)”^’', (129) 


bdl 

2tt 



FIG. 9; Scale dependence of the renormalized disorder b in 
the flat phase for fixed /o = 20 and for <5=1, 5, 25, 100, 300, 
1000, 10000, 100000, increasing from top to bottom. Critical 
curve (<5 = 0) is shown in red. 


A 

bo 



FIG. 10: Scale dependence of the renormalized disorder b 
(measured in the units of its bare value bo) for ycolyccr = 
30 and for different value of bod^/2'K= 0.001, 0.5, 0.9, 1.2, 
1.34, 1.36, 1.37525, 1.382, 1.385, 1.386, 1.38672, 1.3867555, 
1.38676, 1.38677, 1.3868, 1.38685 increasing from bottom to 
top. Critical curve (S — 0, bo = 1.38675444319) is shown in 
red. 


for region (V). Hence, the ripple scale can be estimated 
as 


Lr 



ei/2;. 


for region (IV), 
for region (V). 


(130) 


While r\ is small (r\ = 2/(ic) in the limit of high spatial 
dimensionality of the embedding space, it is a number of 
order unity, rj ~ 0.7 -F 0.8, for a physical 2D membrane 
(e.g., graphene) in a 3D space. This yields the character¬ 
istic scale for ripples of the order of Ginzburg length<2^ 
The functions k{z) and b{z) are plotted in Figs. |8]and 
IHfor /o = 20 and different values of 6. These figures nicely 
illustrate the non-monotonous scale dependencies of the 
coupling constants, with the bending rigidity k(z) having 
a minimum and the disorder b{z) showing a maximum, 
in full agreement with the above analytical results. 
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FIG. 11: Dependence of the parameter 5 (which labels RG 
flow lines) on the bare disorder bodll2'iT for fixed bending 
rigidity, >co/>Cci = 30. The sign change of S corresponds to 
the crumpling transition. 


We also plot the dependence of disorder on the spatial 
scale for a fixed value of >cq = 30>fcr and various values of 
the bare disorder bo (see Fig. uni). In this case, starting 
points of the RG flow lie on a vertical line in Fig. [51 This 
plot serves as a nice illustration of the dependence of the 
characteristic scale Lr on the bare disorder bo. Indeed, 
it is seen, that for small bo the disorder drops quickly at 
the Ginzburg scale, in agreement with Eq. (11301) . With 
increasing bo the length Lr increases. This effect becomes 
particularly strong when disorder becomes comparable 
to the critical one, i.e., the parameter S approaches zero, 
as predicted by Eq. (11271) . For the near-critical curves 
{6 small but still positive) and the critical one (5 = 0), 
the disorder slowly increases up to a parametrically large 
scale Lr and then drops down. On the other side of the 
transition (5 < 0), the disorder shows the same behaviour 
for a while but eventually starts to increase rapidly, which 
reflects the crumpling. In Fig. [11] we plot the parameter 
5 as a function of the bare disorder. 


D. Dynamic and static correlation functions 

To characterize dynamic and static fluctuations in 
the membrane, we introduce the following correlation 
functions)^ 

= Sr,pG^+^{x) (131) 

(/i“(0))(h/3(x)) = 5„/3 G"(x), (132) 

where angular brackets denote the Gibbs averaging, while 
the overline stands for the disorder averaging. The func¬ 
tion , incorporates both dynamical and static fluctu¬ 
ations, while G® includes static correlations only. These 
functions depend on the absolute value of the distance 


X = |x|, so that their Fourier transforms 

G^+® = J 

= 2Tr J G'^'^^{x)Jo{qx)xdx, (133) 

Gq = J G®(a;)e-*'i’^d^x 

= 2Tr J G^{x)Jo{qx)xdx (134) 

depend on the absolute value of momentum q = |q|. 

The dynamic part of the fluctuations is thus given by 
the difference of these two functions, 

= G^+® - G® (135) 

(this function was used in the previous sections without 
index d). The correlation functions defined above can 
be straightforwardly calculated on the basis of the above 
RG analysis. Specifically, we first renormalize the the¬ 
ory from the original ultraviolet scale to the scale 1/q. 
As a result, all non-linear effects get incorporated in the 
renormalization of x and b. Having renormalized the 
couplings, we evaluate the correlation functions at the 
Gaussian level, which yields 


^d-\-s _ 

4 (T + A, 

(136) 


r \>iq J 


Gt = 

T 

(137) 


Xgq'^ 


G% = 

bg 

9^’ 

(138) 


where Xg and bg depend on q according to the RG equa¬ 
tions derived above. 

At this point, it is worth recalling that the running 
scale A for renormalization of x and b is associated with a 
wave-vector q conjugated to the coordinate x in the refer¬ 
ence plane. From the experimental point of view, a more 
natural coordinate on the membrane surface is given by 
the vector r in the embedding space [see Eq. (IT^ ]. As 
seen from the Fig. [T^l a difference between vectors r and 
X becomes stronger when the system approaches critical¬ 
ity. 

The physical correlation functions depend on the dis¬ 
tance r = |r| in the embedding space between two points 
on the membrane surface: 

5'^+®(r) = G^+^{r/0 = |4+Vo(Qr)^,(139) 
9^{r) = G®(r/e) = J ff^Jo(Qr)^, (140) 

/(r) = /+®(r)- 5 ®(r). (141) 

Here ^ is the running stretching parameter that 

relates r and x : r = [see Eq. (IM)) ]. The correlation 














18 


membrane without fluctuations 

L 

-► 


membrane deep in the flat phase 


L 



critical membrane 


R^^lL 



FIG. 12: Top: Membrane with intrinsic size L in the ab¬ 
sence of fluctuations; Center. Membrane deeply in the flat 
phase; Bottom: Membrane at criticality. The linear size, R, 
of membrane with fluctuations is smaller by a factor The 
distance, r, in the embedding space between two points on 
membrane surface scales with the distance on the reference 
plane with a local stretching factor: r ~ Membrane 

shows self-similarity at different scales, so that the dashed line 
represents a coarse-grained shape of the membrane at largest 
scale shown in the picture. The fractal structure of mem¬ 
brane is illustrated in the grey box magnifying a segment of 
the membrane. 


function gg (measured experimentally) is related to Gq 
as follows: 

9QR =i"GQ, (142) 


scaling. Therefore, the difference between g and G is im¬ 
material. On the other hand, for a critical or near-critical 
membrane this difference is of crucial importance. In par¬ 
ticular, the RG flow is controlled by rescaled couplings 
><• and b (rather than by and b.) Further, we find from 
Eq. (11421) that the correlation functions in the embedding 
space, gq, and gg, are given by Eqs. (11361) . (I137L 
and (I138|) with >c and b replaced by >c and 6, respectively. 

In the flat phase, we get the following asymptotic scal¬ 
ing behaviour of the static and dynamic correlation func¬ 
tions at Q —0: 


9q ^ 

1 

(143) 

9q « 

1 

(144) 

Q4-2r, • 


It is worth noting that an analogous asymptotic relation 
between the dynamic and static correlation functions was 
obtained in Ref. for £) = 4 — e for the flat-phase fixed 
point (called P4 point there). Equations (11431) . (11441) 
imply, in particular, that the characteristic dynamic and 
static transverse excursions of a membrane (root- mean- 
square values of the corresponding fluctuations of h) scale 
with its size i? oc L as follows: 


C, cx , (145) 

cx . (146) 

At the crumpling transition point we have = >^cr = 
const, b cx and oc L~'^G ^ Therefore, the char¬ 

acteristic magnitudes of the transverse excursions of the 
membrane scale at the transition as 

oc R, (147) 


and R oc Since >Ccr is proportional to [see 

Eq. (IM)) ] the dimensionless coefficient in Eq. (11471) turns 
out to be small, on the order of 1/ \fdc. 


E. Ripple intensity and correlations 

It is natural to characterize the intensity and spatial 
correlations of ripples (i.e., of static transverse deforma¬ 
tions) with the static dimensionless correlation function 
of spatial gradients of out-of-plane displacements. Such 
function can be expressed in terms of as follows 

R(r) = I ^Q^g^Qj,{Qr). (149) 


with ^ = ^g. In the flat phase, the rescaling factor ^ re- Using Eq. (I142|), we rewrite Eq. (11491) as 
mains finite in the infrared limit (i.e. the membrane lin¬ 
ear size R = in the embedding space is proportional 
to its intrinsic size L), and does not essentially affect the 


H{r) = 


b, 


(27r)2 q2 


9 giqx 


(150) 
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In b 



cay of H{r) changes because of stronger suppression of 
the disorder by increased bending rigidity. Note that at 
very small distances, r < a (where a is the ultraviolet 
cutoff length, which is of the order of the lattice con¬ 
stant), one should replace r with a in the first line of 
Eq. ()151|) . We thus find that the dimensionless parame¬ 
ter H{0) which represents the averaged squared surface 
tilt is given by 


The length 1/q* and the averaged squared tilt H{0) are 
two natural parameters that characterise the character¬ 
istic extension and magnitude of ripples. Importantly, 
both of them decrease with increasing temperature. 


2. Weak disorder [region (V)]. 


FIG. 13: Schematic dependence of the effective disorder 

strength bq on the length scale q~^ on the log-log scale for 
regions (IV) and (V). The indices r;/4 and 2ri denote the 
power-law decay exponents in the corresponding regimes. 


It is worth noting that function H (r) characterizes fluc¬ 
tuations and correlations of normal vectors to the mem¬ 
brane surface. Such fluctuation of tilt angle of the surface 
are directly studied in graphene experiments, see below. 

Let us now discuss the temperature dependence of the 
size and height of the ripples deeply in the flat phase, 
i.e. in regions (IV) and (V) in Fig. [71 In this situation, 
the difference between r and x coordinates is not par¬ 
ticularly important and can be discarded. Using results 
of Secs. IV A[ IV Bl and IV Cl we find the behavior of bq 
(shown schematically in Fig. [T51) and thus for H{r) for 
these two cases. 


1. Strong disorder [region (IV)[. 


The calculation of the integral entering Eq. (I150|) yields 



In Eq. (11511) . 


r < 1/q*, 


1/q* <r<Li^\ 


2r, 


< r. 


(151) 




Q^fol'irt — J_ g46ox:o/3Tr; 

q* q* 


(152) 


is the spatial scale determined by the condition / ~ 1 . 
For r > L^\ the exponent characterising the spatial de¬ 


In the case of a weak disorder, the only difference is 
that disorder falls with the exponent 2irj from the very 
beginning of renormalization. Hence, similar to the case 
of strong potential the characteristic ripple size is deter¬ 
mined by the Ginzburg scale. For D = 2 this scale is 
determined from the condition SA/’Hq ^ 1 [see Eq. (1551) ]. 
yielding 


q =qjg-- 


/3H2Ai(2Ai + 2A)T 
{2fj, + X)yc^ 


(154) 


Therefore, we find the ripple size 

Lj. — — oc (155) 

q* Vt 

The averaged squared surface tilt is given by Eq. (11531) 
which can be rewritten as 

«(») = :!'"(¥)• TO 


where T* ~ h'^{2^x + A)/[ 6 H 2 /r(/i -|- A)a^]. We see that, 
also in this case, both the size of the ripples and the 
characteristic surface tilt a/ 77(0) decrease with increas¬ 
ing temperature. 


3. Comparison with experiment 

Let us compare our results with available experimen¬ 
tal data on suspended graphene. In Ref. [O the param¬ 
eters of ripples at room temperature were found to be 
Lr — 5-F 10 nm for the ripple size and yjHifi) ~ 5° ~ 0.1 
for the characteristic tilt angle. Very similar results 
at r = 300 K were obtained in Ref. HH which found 
Lj. ~ 10 nm and 77(0) ~ 0.01. The small value of 77(0) 
indicates that the system is in the weak-disorder regime. 
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so that the results of Sec. IVE2I are expected to be ap¬ 
plicable. Indeed, the measured room-temperature ripple 
size Lr agrees well with Eqs. (I155|) and (115411 that yield 
Lr ~ 5 nm at T = 300 K. Comparing (11561) with the ex¬ 
perimentally measured H{Q) ~ 0.01, we get an estimate 
for the disorder strength, hy, ~ 0.03. 

The authors of Ref. provided also some informa¬ 
tion about the temperature dependence of the ripple 
characteristics. Specifically, they performed measure¬ 
ments also at T = 150 K and found that, in compari¬ 
son with the room temperature, the ripple size increased, 
Lj. ~ 18 nm whereas the averaged squared surface tilt 
H{Q) remained almost unchanged. These findings are 
in reasonable agreement with our weak-disorder results 
()154|) and ()156|) which predict a square-root increase of 
Lr and a slow (logarithmic) increase of H{Q) with inverse 
temperature. 

It is worth reminding the reader at this point that 
static ripples coexist with dynamical fluctuations. In 
the weak-disorder regime, the relative strength of the 
two types of fluctuations is controlled by the parame¬ 
ter /o = boKo/T, see Eqs. ()137|1 and (1138(1 . According 
to the above estimate, this parameter is close to unity at 
room temperature for the samples experimentally studied 
in Refs. [l3 and 1^. However, it has a 1/T temperature 
dependence which reflects the fact that the dynamical 
fluctuations become suppressed with lowering te mper- 
ature. Contrary to this, the typical tilt angle y^H{0) 
characterising static ripples gets enhanced with lower¬ 
ing temperature, Eq. (11561) . Therefore, measurement of 
the temperature dependence of y^H{0) may be useful for 
experimentally differentiating between the static and dy¬ 
namic fluctuations. 


VI. SUMMARY AND OUTLOOK 

In this article, we have discussed the rippling and the 
crumpling transition in graphene with a static quenched 
disorder. We derived RG equations, Eqs. (l68)l . (|69ll . and 
(EOl) for a model of a crystalline membrane with out- 
of-plane (random curvature) disorder. Equations (1691) 
and (1701) describe a combined flow of the running dimen¬ 
sionless bending rigidity jir/T and the running disorder 
strength b = fTj'k. They yield, in particular, a critical 
curve b{}<) separating the flat and the crumpled phases, 
see Figs. (H [5l and (HI Equation (1551) controls the spa¬ 
tial contraction of the membrane due to its deformation. 
Even deep in the flat phase, random fluctuation of mem¬ 
brane tension caused by the disorder may strongly affect 
the behavior of the bending rigidity ^q/T. Specihcally, 
for a sufficiently strong disorder, the bending rigidity de¬ 
creases at the first stage of the renormalization, reaches 
minimum, and only then starts to grow (Fig. [8]) . Fur¬ 
thermore, we have found that disorder bq also changes 
non-monotonously in the flat phase if the bare disorder 
bo is sufficiently strong. Specifically, bq hrst increases 
slowly (logarithmically) with L, then reaches a maximum 


at a certain scale Lr, and finally decreases according to 
a power law at larger scales, see Figs. iQlandlTOl 

The random static out-of-plane fluctuations of the 
graphene membranes can be identified with experimen¬ 
tally observed ripples, with the length scale Lr play¬ 
ing the role of the characteristic ripple size. The 
found values and temperature dependencies of the ripples 
parameters—the size Lr and the typical surface tilt angle 
—are in a goo d agreement with experimental ob¬ 
servations of Refs. [ITII^ if a disorder strength bo ~ 0.03 
is assumed. 

We have also briefly discussed an in-plane disorder and 
showed that it is irrelevant in the RG sense (if one ex¬ 
cludes a long-range disorder whose correlation function 
is highly singular at small momenta) and thus does not 
affect our main conclusions. The effect of the in-plane 
disorder at atomic scales may, however, be important for 
determining the bare value of the out-of-plane disorder. 

Before closing the paper, we discuss some of possible 
directions of future research. 

1. In our work we considered a free-standing mem¬ 
brane without tension. On the other hand, the 
tension may become essential under certain exper¬ 
imental conditions. In a clean case, such a mem¬ 
brane might demonstrate a buckling transitioni^ 
It remains to be explored what will be the effect of 
disorder in this situation. 

2. We have assumed that the disorder is of short- 
range character. On the other hand, a finite density 
of topological defects (like dislocations or disclina- 
tions) may yield long-range-correlated disorder. An 
earlier workS predicts a variety of possible phases 
in a membrane with long-range disorder. A study 
of crumpling transition and of rippling in a mem¬ 
brane with physically relevant long-range disorder 
remains an interesting prospect for future. 

3. Our analysis did not include terms preventing self¬ 
crossing of membrane which are known to become 
important in the crumpled phase. It remains to 
be investigated whether such terms may affect the 
physics in the near-critical regime. 

4. There is a certain analogy between the physics of 
a membrane and that of Anderson metal-insulator 
transition in disordered (and possibly interacting) 
systems2ii2^. In particular, the field theory of a dis¬ 
ordered membrane developed above bears similar¬ 
ity with the cr-model description of the Anderson 
localization. Within this analogy, the flat phase 
corresponds to a metal, the crumpled phase to an 
insulator, and the dimensionless bending rigidity 
>c/T to the dimensionless conductance. Remark¬ 
ably, both problems manage to evade the Mermin- 
Wagner theorem, showing a transition also in D = 
2. (In the case of Anderson transition, this requires 





















21 


either spin-orbit coupling or electron-electron in¬ 
teraction.) Static fluctuations of local deforma¬ 
tions (ripples) in a disordered membrane can be 
viewed a counterpart of mesoscopic fluctuations of 
wave functions (or local density of states) in the 
Anderson-localization problem. An interesting and 
important question is whether this analogy can be 
pushed further and, in particular, whether the rip¬ 
ple statistics at the crumpling transition is charac¬ 
terised by multifractality that is a hallmark of the 
Anderson-transition critical point. 


the atomic mass M and therefore is large compared to 
typical electronic mass scales. In other words, the l.h.s. 
of Eq. (lAlIl should be larger than its r.h.s. by a factor 
of the order of M/m ^ 2 • 10'’’ (where m is the elec¬ 
tron mass). Hence, for the problem discussed here, the 
flexural phonons can be treated semiclassically for any 
temperature and with a very good precision. The corre¬ 
sponding criterium for longitudinal phonons reads 


T > Tin = 6A2 


+ A) 


(A2) 


5. On the experimental side, a more systematic study 
of rippling and crumpling in free-standing graphene 
would be highly desirable. In particular, measure¬ 
ments of the temperature dependence of ripple pa¬ 
rameters in a broader temperature range would be 
of great interest. Furthermore, experiments on var¬ 
ious kinds of membranes (including emerging 2D 
materials) are expected to be instrumental for ex¬ 
ploring the whole phase diagram of the problem. 


The r.h.s. of this inequality can be estimated as 
(m/M)Ea^ where Ea is a characteristic atomic energy 
scale. Taking Ea = 10 eV, we get a rough estimate 
Tin ~ 5 K. Using known results for the mechanical pa¬ 
rameters of graphene yields a somewhat larger value. 
Tin — 80 K, which still leaves enough room for the valid¬ 
ity of the semiclassical theory. Furthermore, our theory 
remains applicable also at lower temperatures, T < Tin, 
where it describes the physics on sufficiently large spatial 
scales, <7 < T /hs. 
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Appendix A: Applicability of the quasiclassical 
approximation 

In this Appendix, we provide a justification for the 
quasiclassical approximation used in this paper and de¬ 
termine the regime of its validity. 

The quasiclassical approximation is valid for not too 
low temperatures (see also Refs EHHi] for discussion). 
More specifically, it is well justified provided that tem¬ 
perature is large compared to frequencies of both out-of- 
plane and in-plane phonons: T > feuq, T > hsq^ where 
Wq is given by Eq. ([2]), and s = is the 

velocity of the longitudinal in-plane phonons (Here we 
take account that velocity of transverse in-plane phonons 
is smaller for graphene parameters.) The characteristic 
momentum q of the discussed problem is the Ginzburg 
scale g* which itself depends on temperature. According 
to Eq. (II54I) . the condition T > fiwq* can be rewritten as 

Note that this inequality does not contain temperature. 
Substituting graphene parameters (p ~ 7.6 x 10“^ kg/m^, 
A ~ 3 eV/A^, p ~ 9 eV/A^, and k eV), we find that 
l.h.s. of this inequality exceeds the r.h.s. by a factor 
of the order of 10^, so that this requirement is perfectly 
met. This result is not surprising, because the density 
of graphene p, entering this estimate is proportional to 


Appendix B: Screening of h} interaction 

In this Appendix, we present technical details of cal¬ 
culation of the screening. As a starting point we use 
equation for free energy derived in Ref. 

^ ^ Jidk)k‘^\hk\‘^ 4^ / ^dkidk2dkz)Rafijs{q) 

X fciafc2/3A:3T,fc45(hkihk2)(hk3hk4), (BI) 

where q = ki -|- k2 and ki -|- k2 -f- ks -f k4 = 0, and the 
interaction kernel reads 

N 

+ -2- 

Choosing ki = k -f q, k2 = — k, ks = — k' — q, k4 = k', 
we find, after some algebra, 

.^a/37(5 (q)^la^2/3^37^45 — .^q(k, k ), (B3) 

where i?q(k, k') is given by Eq. (l27l) of the main text. 
Screened interaction Rap-yS obeys^ 

Ra0yS — RaPjS RoL^'y'5'^y'5'a'p'Ra.'p'^5-; 
with the tensor polarization operator IV-ySap given by 

= J{dk) kakpk^ksG'lG\_^. (B5) 

Because of the rotation symmetry one can write 

n^,5a^ — H” Bq 

+ {5ap q^qs H-) nW + qaqpqjqs H^^^ (B6) 
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where functions Ilq, Ilq^^ and Ilq^^ depend on |q| only, 

and (H-) stands for sum over permutations of a, (3, 7 , 

and 6. Due to the projection operators P entering in the 
unscreened coupling, Eq. (IB2I) . functions Ilq ^ and Ilq^ 
drop out from Eq. (IB4I) . The function Ilq can be easily 
obtained by multiplication of Eq. (jB5p by PapP-yS and 
summation over repeated indices. Taking into account 
that trace of the matrix Pap equals to D — 1 , we arrive at 
Eq. (1571) of the main text. Further, substituting Eq. (IB5I) 
into Eq. (IB4p , we find that the screened interaction Rap-yS 
can be written in the same form as Eq. (IB2I) but with Af 
and M replaced by their screened values A/’q and A4q. 
These values are given by Eqs. (l33l) and (l34l) of the main 
text, respectively. 


Appendix C: Regularization integrals 

In order to shed light on a connection of our ap¬ 
proach to the SCSA, we employ the regularization in¬ 
tegrals A{D, rj) and B{D, rj) introduced in Ref. IH: 

/(‘‘" b'-lx I nr-- 

22D+V+1 yy{D-l)/2 p2 ^4^^ I 

BiD,,) = (C2) 

_ _ i)r r (f) r(2 - 77) 

22+d ^d/2 y ( 4 ^) r r ’ 



Here n is an arbitrary unit-length vector and xj^ = x — 
n(xn). These integrals naturally arise when one uses for 
calculation of the polarization loop the Green function, 
Eq. (H71) . with the self-energy found self-consistently by 
replacement of the bare Green function in Eq. (l43l) 
with Gk. Evaluation of both integrals can be performed 
in the following way. First, one uses the identity 


1 1 
x^lx- “ r(|)r(^)' 

00 CXD 

X J dh J dt2 

0 0 


where 0 = 4 — 77 for the integral A and 9 = D + 2ri 
for the integral B. The integral over dx = dx^/( 27 r)'° 
becomes then a Gaussian one and is easily calculated. 
The remaining integral can be done by using the following 
change of variables: ty = re^ coshz and ±2 = re“^ coshz. 
After lengthy but straightforward calculations one arrives 
at Eqs. (ICll) and (IG2D . 


^ K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. 
Zhang, S.V. Dubonos, I.V. Grigorieva and A.A. Firsov, 
Science 306, 666 (2004). 

^ K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, 
M.I. Katsnelson, I.V. Grigorieva, S.V. Dubonos and A.A. 
Firsov, Nature 438, 197 (2005). 

^ Y. Zhang, Y.-W. Tan, H.L. Stormer and P. Kim, Nature 
438, 201 (2005). 

4 A.K. Geim and K.S. Novoselov, Nature Materials 6, 183 
(2007). 

® A.H. Gastro Neto, F. Guinea, N.M.R. Peres, K.S. Novo¬ 
selov, and A.K. Geim, Rev. Mod. Phys. 81, 109 (2009). 

® S. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, 
Rev. Mod. Phys. 83, 407 (2011) 

^ V.N. Kotov, B. Uchoa, V.M. Pereira, F. Guinea, and 
A. H. Castro Neto, Rev. Mod. Phys. 84, 1067 (2012) 

® M.I. Katsnelson, Graphene: Carbon in Two Dimensions 
Hardcover, Cambridge University Press (2012) 

® E. L. Wolf, Graphene: A New Paradigm in Condensed 
Matter and Device Physics, Oxford University Press (2014) 
4° L.E.F. Foa Torres, S. Roche, J.-C. Charlier, Introduction 
to Graphene-Based Nanomaterials From Electronic Struc¬ 


ture to Quantum Transport, Cambridge University Press 
(2014) 

44 K.S. Novoselov, Z. Jiang, Y. Zhang, S.V. Morozov, H.L. 
Stormer, U. Zeitler, J.C. Maan, G.S. Boebinger, P. Kim, 
and A.K. Geim, Science 315, 1379 (2007). 

4^ K.I. Bolotin, K.J. Sikes, Z. Jiang, M. Klima, G. Fudenberg, 
J. Hone, P. Kim, and H.L. Stormer, Solid State Commun. 
146, 351 (2008). 

4® X. Du, I. Skachko, A. Barker, and E.Y. Audrey, Nature 
Nanotechnology 3, 491 (2008). 

44 K.I. Bolotin, K.J. Sikes, J. Hone, H.L. Stormer, and P. 

Kim, Phys. Rev. Lett. 101, 096802 (2008). 

4® X. Du, 1. Skachko, F. Duerr, A. Luican, and E.Y. 

Andrei, Nature 462, 192 (2009). 

4® K.I. Bolotin, F. Ghahari, M.D. Shulman, H.L. Stormer, 
and P. Kim, Nature 462, 196 (2009). 

44' J.C. Meyer, A.K. Geim, M.I. Katsnelson, K.S. Novoselov, 
T.J. Booth and S. Roth, Nature 446, 60 (2007). 

4® J. Scott Bunch, A.M. van der Zande, S.S. Verbridge, I.W. 
Frank, D.M. Tanenbaum, J.M. Parpia, H.G. Craighead, 
and P.L. McEuen, Science 315, 490 (2007). 

4® F. Miao, S. Wijeratne, Y. Zhang, U.C. Coskun, W. Bao, 














23 


and C.N. Lau, Science 317, 1530 (2007). 

R. Danneau, F. Wu, M.F. Craciun, S. Russo, M.Y. Tomi, 

J. Salmilehto, A.F. Morpurgo, and P.J. Hakonen, Phys. 
Rev. Lett. 100, 196802 (2008). 

C. Gomez-Navarro, M. Burghard, and K. Kern, Nano Lett. 
8, 2045 (2008). 

D. C. Elias, R.V. Gorbachev, A.S. Mayorov, S.V. Mo¬ 
rozov, A.A. Zhukov, P. Blake, L.A. Ponomarenko, I.V. 
Grigorieva, K.S. Novoselov, F. Guinea, and A.K. Geim, 
Nature Physics, 7, 701 (2011). 

2^ C.N. Lau, W. Bao, and J. Velasco, Materials Today 15, 
238 (2012). 

2^ D. Nelson, T. Piran, and S. Weinberg (Eds.) Statistical 
Mechanics of Membranes and Surfaces (World Scientific, 
Singapore, 1989). 

2® L.M. Woods and G.D. Mahan, Phys. Rev. B 61, 10651 

( 2000 ). 

2® H. Suzuura and T. Ando, Phys. Rev. B 65, 235412 (2002). 
2^ E.H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 
(2007). 

2® J.L. Manes, Phys. Rev. B 76, 045430 (2007). 

2® A.H. Castro Neto and E.A. Kim, Euro. Phys. Lett. 84, 
57007 (2008). 

A. Fasolino, J.H. Los, and M.I. Katsnelson, Nature Mate¬ 
rials 6, 858 (2007). 

D. M. Basko and I.L. Aleiner, Phys. Rev. B 77, 041409(R) 
(2008). 

®2 D.M. Basko, Phys. Rev. B 78, 125418 (2008). 

E. Mariani and F. von Oppen, Phys. Rev. Lett. 100, 
076801 (2008). 

F. von Oppen, F. Guinea, and E. Mariani, Phys. Rev. B 
80, 075420 (2009). 

E. Mariani and F. von Oppen, Phys. Rev. B 82, 195403 

( 2010 ). 

M.A.H. Vozmediano, M.I. Katsnelson, and F. Guinea, 
Physics Reports 496 109, (2010). 

E.V. Castro, H. Ochoa, M.I. Katsnelson, R.V. Gorbachev, 
D.C. Elias, K.S. Novoselov, A.K. Geim, and F. Guinea, 
Phys. Rev. Lett. 105, 266601 (2010). 

K. V. Zakharchenko, R. Roldan, A. Fasolino, and M.I. Kat¬ 
snelson, Phys. Rev. B 82, 125435 (2010). 

R. Roldan, A. Fasolino, K.V. Zakharchenko, and M.I. Kat¬ 
snelson, Phys. Rev. B 83, 174104 (2011). 

P. San-Jose, J. Gonzalez, and F. Guinea, Phys. Rev. Lett. 
106, 045502 (2011). 

H. Ochoa, E.V. Castro, M.I. Katsnelson, and F. Guinea, 
Phys. Rev. B 83, 235416 (2011). 

■^2 I.V. Gornyi, V.Yu. Kachorovskii, and A.D. Mirlin, Phys. 
Rev. B 86, 165413 (2012). 

S. Das Sarma and E. H. Hwang, Phys. Rev. B 87, 035415 
(2013). 

M. V. Medvedyeva and Ya. M. Blanter, Phys. Rev. B 88, 
125423 (2013). 

K. S. Tikhonov, Wei L.Z. Zhao, and A. M. FinkeLstein, 
Phys. Rev. Lett. 113, 076601 (2014). 

A. A. Balandin, S. Ghosh, W. Bao, 1. Calizo, D. Tewelde- 
brhan, F. Miao, and C. N. Lau, Nano Lett. 8, 902 (2008); 
A. A. Balandin, Nature Materials 10, 569 (2011). 

H. Seol, 1. Jo, A. L. Moore, L. Lindsay, Z. H. Aitken, M. 

T. Pettes, X. Li, Z. Yao, R. Huang, D. Broido, N. Mingo, 
R. S. Ruoff, and L. Shi, Science 328, 213 (2010). 

Z. Wang, R. Xie, C. T. Bui, D. Liu, X. Ni, B. Li, and J. 
T. L. Thong, Nano Lett. 11, 113 (2011). 

K.H. Michel, S. Costamagna, and F.M. Peelers, Phys. Rev. 


B 91, 134302 (2015). 

N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 

L. D. Landau and E.M. Lifshitz, Statistical Physics, Part I 
(Pergamon Press, Oxford, 1980). 

D. R. Nelson and L. Peliti, J. Phys. (Paris) 48, 1085 (1987). 
Y.Kantor and D.R. Nelson, Phys. Rev. Lett. 58, 2774 
(1987); Phys. Rev. A 36, 4020 (1987); 

M. Paczuski, M. Kardar, and D.R. Nelson, Phys. Rev. Lett. 
60, 2638 (1988). 

F. David and E. Guitter, Europhys. Lett. 5, 709 (1988). 

E. Guitter, F. David, S. Leibler, and L.Peliti, Phys. Rev. 
Lett. 61, 2949 (1988). 

J.A. Aronovitz and T.C. Lubensky, Phys. Rev. Lett. 60, 
2634 (1988). 

E. Guitter, F. David, S. Leibler, and L. Peliti, J.Phys. 
France 50 1787 (1989). 

J. Aronovitz, L. Golubovic, and T.C. Lubensky, J.Phys. 
France 50 609 (1989). 

M. Paczuski and M. Kardar, Phys. Rev. A 39, 6086 (1989) 

L. Radzihovsky and D.R. Nelson, Phys. Rev. A 44, 3525 
(1991). 

D.R. Nelson and L. Radzihovsky, Europhys. Lett. 16, 79 
(1991). 

G. Gompper and D.M. Kroll, Europhys. Lett. 15, 783 
(1991). 

L. Radzihovsky and P. Le Doussal, J.Phys. I France 2 599 
(1992). 

P. Le Doussal and L. Radzihovsky, Phys. Rev. Lett 69, 
1209 (1992). 

D.C. Morse, T.C. Lubensky, and G.S. Grest, Phys. Rev. A 
45, R2151 (1992). 

P. Le Doussal and L. Radzihovsky, Phys. Rev. B 48, 3548 
(1993). 

M. J. Bowick, S.M. Catterall, M. Falcioni, G. Thorleifs- 
son, and K.N. Anagnostopoulos, J. Phys. I France 6, 1321 
(1996). 

J.-P. Kownacki, and D. Mouhanna, Phys. Rev. E 79, 
040101(R) (2009). 

D. Gazit, Phys. Rev. E 80, 041117 (2009). 

D. Gazit, Phys. Rev. B 80, 161406(R) (2009). 

F. L. Braghin and N. Hasselmann , Phys. Rev. B 82, 
035407 (2010). 

V.V. Lebedev and E.I. Kats, Phys. Rev. B 85, 045416 

( 2012 ). 

E. I. Kats and V.V. Lebedev, Phys. Rev. B 89, 125433 
(2014). 

B. Amorim, R. Roldan, E. Cappelluti, A. Fasolino, F. 
Guinea, and M. 1. Katsnelson, Phys. Rev. B 89, 224307 
(2014). 

E. 1. Kats and V. V. Lebedev, Phys. Rev. E 91, 032415 
(2015). 

A. L. Vazquez de Parga, F. Calleja, B. Borca, M. C. G. 
Passeggi, Jr., J. J. Hinarejos, F. Guinea, and R. Miranda 
Phys. Rev. Lett. 100, 056807 (2008). 

V. W. Brar, Y. Zhang, Y. Yayon, T. Ohta, J. L. McCh- 
esney, A. Bostwick, E. Rotenberg, K. Horn, and M. F. 
Crommie, Appl. Phys. Lett. 91, 122102 (2007). 

V. Geringer, M. Liebmann, T. Echtermeyer, S. Runte, M. 
Schmidt, R. Ruckamp, M. C. Lemme, and M. Morgenstern, 
Phys. Rev. Lett. 102, 076102 (2009). 

W. Bao, F. Miao, Z. Chen, H. Zhang, W. Jang, C. Dames, 
and C. N. Lau, Nature Nanotechnology 4, 562 (2009) 

D. A. Kirilenko, A. T. Dideykin, and G. Van Tendeloo, 


50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 


24 


Phys. Rev. B 84, 235417 (2011) 

L. Tapaszto, T. Dumitrica, Sung Jin Kim, P. Nemes-Incze, 
Chanyong Hwang, and L. P. Biro, Nature Physics 8, 739 
( 2012 ). 

D.A. Kirilenko, Technical Physics Letters 39, 325 (2013). 
P. Xu, M. Neek-Amal, S. D. Barber, J. K. Schoelz, M.L. 
Ackerman, P. M. Thibado, A. Sadeghi, and F.M. Peelers, 
Nature Communications 5, 3720 (2014); M. Neek-Amal, 
P. Xu, J.K. Schoelz, M.L. Ackerman, S.D. Barber, P.M. 
Thibado, A. Sadeghi, and F.M. Peelers, ibid 5, 4962 
(2014). 

J. K. Schoelz, P. Xu, V. Meunier, P. Kumar, M. Neek- 
Amal, P. M. Thibado, and F. M. Peelers, Phys. Rev. B 
91, 045413 (2015). 

C. F. Schmidt, K. Svoboda, N. Lei, 1. B. Petsche, L. E. 
Berman, C. R. Safinya, and G. S. Grest, Science 259, 952 
(1993). 

T. Hwa, E. Kokufuta, T. Tanaka, Phys. Rev. A 44, R2235 
(1991); X. Wen, C.W. Garland, T. Hwa, M. Kardar, E. 
Kokufuta, Y. Li, M. Orkisz, and T. Tanaka, Nature 355, 
426 (1992); G. Cheng and D. Li, Adv. Mater. 25, 13 (2013). 
S. Costamagna, M. Neek-Amal, J. H. Los, and F. M. 
Peelers, Phys. Rev. B 86, 041408(R) (2012). 

R.R. Chianelli, E.B. Prestridge, T.A. Pecoraro, and J.P. 
Deneufville, Science 203, 1105 (1979); J. Brivio, D.T. L. 
Alexander, and A. Kis, Nano Lett. 11, 5148 (2011). 

D. Pacile, J. C. Meyer, C. O. Girit, and A. Zettl, Appl. 


Phys. Lett. 92, 133107 (2008). 

The stability of the SGSA with respect to higher-order 
contributions was addressed in Ref. where the inclu¬ 
sion of second-order corrections in the parameter 1/dc = 1 
was found to yield rj « 0.789 for physical membranes. The 
SCSA has been applied to graphene in Ref. [3^ . 

As follows from Eqs. (O and uni), renormalization of ^ 
starts from the inverse ultraviolet cutoff length, l/o. Here 
we neglect renormalization of ^ in the interval q* < q < 1/a 
assuming that {1/dc) ln(l/q'*a) <C 1. 

The Ginzburg length as a characteristic scale for thermally 
excited ripples was also found in Ref, [t^. Numerical simu¬ 
lations of graphene flakes performed in Ref. also yielded 
the s pat ial scale of ripples close to the Ginzburg length. In 
Refs. and [z3, the size of ripples was related to the cou¬ 
pling of phonons with electrons, yielding a T-independent 
length (ultraviolet atomic scale). Interestingly, taking into 
account the numerical factors, this scale was estimated^ 
to be close to the thermal Ginzburg scale corresponding to 
the room temperature. We emphasize that the membrane 
fluctuations discussed in the above works are dynamical 
rather than static, see Sec. EDI 

F. Evers and A.D. Mirlin, Rev. Mod. Phys. 80, 1355 
(2008). 

A.M. Finkelstein, Int. J. Mod. Phys. B 24, 1855 (2010). 




